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Abstract 

Real-world road networks have an approximate scale-invariance 
property; can one devise mathematical models of random networks 
whose distributions are exactly invariant under Euclidean scaling? This 
requires working in the continuum plane. We introduce an axiomati- 
zation of a class of processes we call scale-invariant random spatial 
networks, whose primitives are routes between each pair of points in 
the plane. We prove that one concrete model, based on minimum- 
time routes in a binary hierarchy of roads with different speed limits, 
satisfies the axioms, and note informally that two other constructions 
(based on Poisson line processes and on dynamic proximity graphs) 
are expected also to satisfy the axioms. We initiate study of structure 
theory and summary statistics for general processes in this class. 
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1 Introduction 

Familiar web sites such as Google maps provide road maps on adjustable 
scale (zoom in or out) and a suggested route between any two specified 
addresses. Given k addresses in a country, one could find the route for each 
of the (2) pairs, and call the union of these routes the subnetwork (of the 
country's entire road network) spanning the k points. 

We abstract this idea by considering, for each pair of points (z, z') in the 
plane, a random route TZ{z,z') = TZ{z' ,z) between z and z' . The collection 
of all routes (as z and z' vary) defines what one might call a continuum 
random spatial network, an idea we explain informally in this introduction 
(precise definitions will be given in section [2^ . 
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In particular, for each finite set (zi, . . . , Zk) of points we get a random 
network span(zi, . . . , z^), the spanning subnetwork hnking the points, con- 
sisting of the union of the routes TZ{zi,Zj). Mathematically natural struc- 
tural properties we will impose on the distribution of such a process are 

(i) translation and rotation invariance 

(ii) scale-invariance. 

For < c < oo the scaling map dc : — )• M? takes z to cz; we emphasize 
that (ii) means "naive Euclidean scaling", i.e. invariance under the action 
of (Tc, not any notion of "scaling exponent". For instance, scale-invariance 
implies that the route-length Dr between points at (Euclidean) distance r 

apart must scale as Dr = rDi, where of course 1 < Di < oo. The setup 
so far does not exclude the possibility that routes are fractal, with infinite 
length, and such cases do in fact arise naturally in the tree-like models of 
section [5221 But, envisaging road networks rather than some other physical 
structure, we restrict attention to the case EDi < oo. There is a rather 
trivial example, the complete network in which each TZ{zi, Z2) is the straight 
line segment from zi to Z2, but the assumption "^ < 00" below will exclude 
this example. 

Much of our study involves sampled spanning subnetworks 5(A), as fol- 
lows. Write H(A) for a Poisson point process in of intensity A, independent 
of the network. Then the points ^ of H(A), together with the routes TZ{C, C') 
for each pair of such points, form a random subnetwork we denote by 5(A). 
The distribution of 5(A) inherits the properties of translation- and rotation- 
invariance, and a form of scale-invariance described at ()12p . In particular 
we can define a constant < ^ < 00 by 

i = mean length-per- unit-area of 5(1) 

(where "mean length-per- unit-area " is formalized by edge-intensity at ([T])). 
In section [531 we note a crude lower bound £> \. We impose the property 

^ < 00. 

Regard i as "normalized network length", for the purpose of comparing 
different networks. 

Everything mentioned so far makes sense when only finite-dimensional 
distributions TZ{zi,Zj) are specified. A first context in which we want to 
consider a process over the whole continuum concerns the following conve- 
nient abstraction of the notion of "major road". Write 7^(i)(2;i, 2:2) for the 
part of the route Tl{zi,Z2) that is at distance > 1 from each of zi and Z2- 
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Conceptually, we want to study an edge-process E viewed as the union of 
7^(x)(zi,Z2) over all pairs {z\,Z'i). To formalize this directly would require 
some notion of "regularity" for a realization, for instance some notion of 
a.e. continuity of routes TtiyZx^Z'i) as z\ and z^ vary. But we can avoid this 
complication by first considering only z\,Z2 in H(A) and then letting A — )• (X). 
After defining E in this way, we can define 

:= mean length-per-unit-area of E 

and impose the requirement 

< oo. 

If a process of random routes z') satisfies the properties we have de- 
scribed (as stated precisely in section I2.2p . then we will call it a scale- 
invariant random spatial network (SIRSN). As the choice of name suggests, 
it is the scale- invariance that makes such processes of mathematical interest; 
in section [LSl we briefiy discuss its plausibility for real- world networks. 

We do not know any closely related previous work. We will discuss one 
related area of theory (discrete random spatial networks; section II. 2p and 
one area of application (fast algorithms for shortest routes:; section II. 4p . 
Several more distantly related topics are mentioned in section [821 

1.1 Outline of paper 

The purpose of this paper is to initiate study of SIRSNs, with three em- 
phases. First, we give a careful formulation of an axiomatic setup for 
SIRSNs, with discussion of possible alternatives (section [2]). Second, it 
is not obvious that SIRSNs exist at all! We give details of one construc- 
tion in section [3l That construction envisages a square lattice of freeways, 
with "speed level j" freeways spaced 2-' apart, and the routes are the mini- 
mum time paths. Being based on the discrete lattice makes some estimates 
technically straightforward, but completing the details of proof requires sur- 
prisingly intricate arguments. This construction is somewhat artificial in not 
naturally having all the desired invariance properties, so these need to be 
forced by external randomization. We briefly mention two other construc- 
tions (in section 14.11 based on a weighted Poisson line process representing 
the different-level freeways, and in section |4?2] based on a dynamic construc- 
tion of random points and roads added accoding to a deterministic rule) 
which intuitively seem more natural but for which we have been unable to 
complete all the details of a proof. 
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Third, in sections [5] - [6] we begin developing some general theory from 
the axiomatic setup. Of course scale-invariance is a rather weak assump- 
tion, loosely analogous to stationarity for a stochastic process, so one cannot 
expect sharp results holding throughout this general class of process. Our 
general results might be termed "structure theory" and concern existence 
and uniqueness issues for singly- and doubly-infinite geodesies, continuity 
of routes Tl{zi,Z2) as a function of {zi,Z2), numbers of routes connecting 
disjoint subsets, and bounds on the parameters KDi,£,p{l). One feature 
worth emphasis is that (very loosely analogous to entropy rate for a station- 
ary process) the quantity p{l) is a non-obvious statistic of a SIRSN, but 
turns out to be central in the foundational setup, in the structure theory, 
and in conceptual interpretation as a model for road networks. The latter 
is best illustrated by the "algorithms" story in sections 11.41 and 16.41 

Being a new topic there are numerous open problems, both conceptual 
and technical, stated in a final discussion section [8l 

Before starting technical material, sections 11.21 - 11.51 give further verbal 
discussion of background to the topic. 

1.2 Discrete spatial networks 

Traditional models of (deterministic or random) spatial networks start with 
a discrete set of points and then assign linking edges via some rule, e.g. the 
random geometric graph [20j or proximity graphs [15], surveyed in [5]. One 
specific motivation for the present work was as a second attempt to resolve 
a paradox - more accurately, an unwelcome feature of a naive model - in 
the discrete setting, observed in [8] . In studying the trade-off between total 
network length and the effectiveness of a network in providing short routes 
between discrete cities, one's first thought might be to measure the latter 
by the average, over all pairs {x,y), of the ratio 

(route-length for x to y)/ (Euclidean distance from x to y) 

instead of averaging over pairs at Euclidean distance r to get our KD^. But 
it turns out that (in the n — )• oo limit of a network on n points) one can make 
this ratio tend to 1 for a network whose length is only l-|-o(l) times the length 
of the Steiner tree, by simply superimposing on the Steiner tree a sparse 
Poisson line process. Such "theoretically optimal" networks are completely 
unrealistic, so there must be something wrong with the optimization criteria. 
What's wrong is that the networks are ineffective for small r. One way to 
get a non-trivial tradeoff in the n — )• oo limit was described in [5]: using 
the statistic max^ r~^E,Dr as the measure of effectiveness leads to more 
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realistic-looking networks. In the discrete setting a network model cannot be 
precisely scale- invariant, but such considerations prompted investigation of 
continuum models which are assumed to be scale- invariant, so that r~^E,Dr 
is constant. 

We emphasize that our networks involve roads at definite positions in 
the plane. There is substantial recent literature, discussed in [10], involv- 
ing quite different notions of random planar networks, based on identifying 
topologically equivalent networks. 

1.3 Visualizing spanning subnetworks 

Both construction and analysis of general SIRSNs are based on studying 
subnetworks span(2;i, . . . ,Zk) on fixed or (most often) random points. It is 
helpful to visualize what subnetworks look like - see Figure 1. 



Figure 1. Schematic for the subnetwork of a SIRSN on 7 points • 

The qualitative appearance of Figure 1 is quite different from that of familiar 
spatial networks mentioned above, based on a discrete set of points, which 
could be viewed as abstractions of an inter-city road network, with cities as 
points. In contrast, we are abstracting the idea of the points • being indi- 
vidual street addresses a long way apart. The real- world route between two 
such street addresses will typically consist, in the middle, of roughly straight 
freeway segments but, nearing an endpoint, of a more jagged trajectory of 
shorter segments of lower-capacity roads; our setup and Proposition 1131 im- 
ply the same behavior in our model. 
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1.4 Very fast shortest path algorithms 

There is an interesting connection with the "shortest path algorithms" ht- 
erature. Onhne mapping services and GPS devices require very quick com- 
putations of shortest routes. In this context, the U.S road network is rep- 
resented as a graph on about 15 milhon street intersections (vertices) with 
edges (road segments) marked by distance (or typical driving time), and a 
given street address is recognized as being between two specific street in- 
tersections. Given a pair of (starting and destination) points, one wants 
to compute the shortest route. Neither of the two extremes - pre-compute 
and store the routes for all possible pairs; or use a classical Dijkstra-style 
algorithm for a given pair without any preprocessing - is practical. Bast 
et al (see [9] for an outline) find a set of about 10,000 intersections (which 
they call transit nodes) with the property that, unless the start and destina- 
tion points are close, the shortest route goes via some transit node near the 
start and some transit node near the destination. Given such a set, one can 
pre-compute shortest routes and route-lengths between each pair of transit 
nodes; then answer a query by using the classical algorithm to calculate 
the route lengths from starting (and from destination) point to each nearby 
transit node, and finally minimizing over pairs of such transit nodes. 

This idea is actually used commercially (and patented). Mathematical 
discussion was initiated by Abraham et al [T], who introduced the notion of 
highway dimension, defined as the smallest integer h such that for every r 
and every ball of radius 4r, there exists a set of h vertices such that every 
shortest route of length > r within the ball passes through some vertex in 
the set. They discuss several algorithms whose performance can be analyzed 
in terms of highway dimension, and devise a particular model (a dynamic 
spanner construction on vertices given by an adversary) designed to have 
bounded highway dimension. 

Now saying one can find h independent of r is a form of approximate 
scale-invariance, so the empirical fact that one can find transit nodes in 
the real-world road networks is a weak form of empirical scale-invariance. 
Within our model where precise scale-invariance is assumed, we can derive 
quantitative estimates relating to transit nodes - see section 16.41 

Incidently, the way we define edge-processes £ = £{X, r) in terms of 
routes (mentioned in the Introduction and defined in section 12. 2p is closely 
related to the notion of reach in the algorithmic literature [14j . 
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1.5 Visualizing scale-invariance 



Visualizing a photo of a road, scale-invariance seems implausible, because it 
implies existence of roads of arbitrarily large and arbitrarily small "sizes" , 
however one interprets "size". But scale-invariance is not referring to the 
physical roads but to the process of "shortest routes", as in the discussion 
above. Figure 2 illustrates one aspect of scale-invariance. There is some 
number of crossing places (over the line) used by routes from one square 
to the other square. In our model, scale-invariance implies that the mean 
number of such crossings does not depend on the scale of the map. One 
could test this as a prediction about real-world road networks. 



Figure 2. Schematic for long-distance routes. 

As another empirical aspect of scale-invariance, [16] studied proportions of 
route-length, within distance-r routes, spent on the i'th longest road seg- 
ment in the route (identifying roads by their highway number designation) 
and observe that in the U.S. the averages of these ordered proportions are 
around (0.40,0.20,0.13,0.08,0.05) as r varies over a range of medium to 
large distances. Again, in our models (identifying roads as straight seg- 
ments) scale-invariance implies there is some vector of expected proportions 
that is precisely independent of r. 

2 Technical setup 

In formulating an axiomatic setup there are several alternative choices one 
could make. In section 12.21 we state concisely the choices we made; section 
12.31 discusses alternatives, reasons for choices, and immediate consequences 
or non-consequences of the setup. 

2.1 Stochastic geometry background 

We quote a fundamental identity from stochastic geometry (see \22\ Chapter 
8). Let £ be an edge process - for our purposes, a union of line segments - 
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whose distribution is invariant under translation and rotation. Then E has 
an edge-intensity, a constant l = intensity(£^) G [0, oo] such that 

£;(length oi S D A) = l x area(A), A c M^. (1) 

Moreover the positions and angles at which £ intersects the x-axis (and 
hence any other line) are such that 

mean number intersections per unit length = 27r~^ x intensity(iS) (2) 

and the random angle 6 G (0, vr) of a typical intersection has density 

fe{e) = ^sme. (3) 

2.2 Definitions 

Here we organize the setup via four aspects. 

Some notation. denotes the origin; disc(2;, r) and circle(2;, r) denote the 
closed disc and the circle centered at z. 

Aspect 1. Allowed routes and route-compatability. Define a jagged 
route between two points z, z' of to consist of straight line segments 
between successive points (zj,— oo < i < oo) with limj_!._oo -^j = z and 
liiRi^oo Zi = z' , and such that the total length \zi — Zj-i] is finite. 

A feasible route is either a jagged route or the variant with a finite or semi- 
infinite set of successive line segments; we further require that the route be 
non-sclf-intcrsccting. Write r(z, z') for a feasible route, which from now on 
we will just call route. Wc envisage a route r{z, z') as a one-dimensional 
subset of M?, equipped with a label indicating it is the route from z to z'. 
The route r{z', z) is always the reversal of r{z, z'). 

When we have a collection of routes, we require the following pairwise 
compatability property. 

If two routes r(zi, Zj), r{z[, Zg) meet at two points then the routes 

coincide on the subroute between the two meeting points. (4) 

Aspect 2. Subnetworks on locally finite configurations. Given a 

locally finite configuration of points (zj) in the plane, and routes ir{zi,Zj) 
satisfying the pairwise compatability property, write s for the union of all 
these routes. If s has the "finite length in bounded regions" property 

len(s n disc(0,r)) < oo for each r < oo (5) 
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then call s a feasible subnetwork. Here "len' denotes "length". Formally s 
consists of the vertex set (zj), an edge set which is the union of the edge sets 
comprising each r(zj, zj), and marks on edges to indicate which routes they 
are in. Inclusion s(l) C s(2) means that s(2) can be obtained from s(l) by 
adding extra vertices and associated routes. 

As outlined in section 12.31 there is a natural cj-field that makes the set of 
all feasible subnetworks into a measurable space, so it makes sense below to 
talk about random feasible subnetworks. 

Aspect 3. Desired distributional properties of subnetv^orks. The 

precise definition of the class of processes we shall study uses "finite-dimensional 
distributions" (FDDs), as follows. Given a finite set zi,...,Zk let fizi,...,Zk 



be the distribution of a random feasible subnetwork span(zi, . . . , Zfc) on 
zi, . . . ,Zk- Suppose a family (indexed by all finite sets) of FDDs satisfies 

the natural consistency condition (6) 
invariance under translation and rotation (7) 
invariance under scaling. (8) 



To be precise about ([8]), recall that the scaling map cjc : — )• M? takes z to 
cz. Then the action of ac on span(2:i, . . . , z^) gives a random subnetwork 
whose distribution equals the distribution of span(cjc2i, . . . ,acZk)- 

Appealing to the Kolmogorov extension theorem, we can associate with 
such a family a process of routes TZ{zi, Z2), for each pair zi, Z2 in M^, though 
for a process defined in that way we can only discuss properties determined 
by FDDs. 

As mentioned earlier, much of our study involves sampled spanning sub- 
networks, as follows. For each < A < cxd let H(A) be a Poisson point 
process of intensity A (we sometimes call this point-intensity to distinguish 
from edge-intensity at ([T|)). Make a process (H(A), < A < 00) by coupling 
in the natural way (take a space-time Poisson point process and let H(A) be 
the positions of points arriving during time [0, A]). Taking Poisson points 
independent of the process of routes, we can define 5(A) as the subnetwork 
of routes TZi^^i') for pairs in H(A). We want the resulting processes 
5(A) to have the following properties. 

for each A, 5(A) is a random feasible subnetwork on vertex-set S(A) (9) 

for each A, 5(A) has translation- and rotation-invariant distribution (10) 

S{\i) C S{\2) for Ai < A2 (11) 

applying CTc to 5(A) gives a network distributed as 5(c^^A). (12) 
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For (fT2]) . recall that applying cJc to H(A) gives a point process distributed 
as S(c~^A). 

We omit full measure-theoretic details of the construction of 5(A), and 
just point out what extra conditions are needed to obtain properties ([9] - 
\T2\i . First, we need to impose the technical condition 

the map {zi, . . . , z^) — s- fJ-zi,...,Zk is measurable (13) 

to ensure that 5(A) is measurable. Second, part of the "feasible" assertion in 
@ is the "finite length in bounded regions" property and this property 
for 5(A) cannot be a consequence of assumptions on FDDs only, so we need 

len(5(A) n disc(0,r)) < oo a.s. for each r < oo (14) 

and this will follow from the stronger assumption (jl6p below. 

Aspect 4. Final definition of a SIRSN. To summarize the above: 
given a process of routes lZ{zi,Z2) with FDDs satisfying iP- [8l [T3l). we can 
define the process of sampled subnetworks (5(A), < A < oo) which, if (fH|) 
holds, will have properties (f9l- I12p . Finally, we define a SIRSN as a process 
(denoted by the routes TZ{zi, Z2) or by the sampled subnetworks (5(A), < 
A < 00)) satisfying these assumptions ([6] - El [T3l) and also satisfying the 
extra conditions (|15l20p below. These extra conditions merely repeat and 
formalize the requirements, stated in the introduction, that certain statistics 
be finite. As noted above, these assumptions imply that ([9l- [T2]) hold. 

Write 1 = (1,0) and Di := len 7^(0,1). So Di represents route-length 
between points at distance 1 apart. Our definition of feasible route implies 
1 < Di < 00 a.s., and we impose the requirement 

1 < EDi < 00. (15) 

Next, our definition oi feasible subnetwork implies that 5(1) must have a.s. 
finite length in a bounded region. We impose the stronger requirement of 
finite expected length. In terms of the edge-intensity ([I]), we require 

£ := intensity (5(1)) < 00. (16) 

Finally, we define 

£{X,r):= [j 7^(C,0\(disc(e,r)Udisc(e^r)) (17) 
€,€'eH(A) 
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and edge-intensities 



p{X,r) 
p{r) 



intensity(iS(A, r)) 
lim p(A, r) 



A— >CXD 



(18) 
(19) 



and impose the requirement 



p{l) < oo 



(20) 



whose significance is discussed in the next section. Lemma[20]will show that 
(pOj) implies (fT6|) . If we do not require ([20]) but instead require (fT6]) , call the 
process a weak SIRSN. 

2.3 Discussion of technical setup 

Aspect 1. Allowed routes and route-compatability. Because we 
want routes to have a well-defined lengths, a minimum assumption would 
be that routes are rectifiable curves. We have assumed "feasible routes" in 
order to simplify notation. We believe that the theory would be essentially 
unchanged if instead one allowed rectifiable curves, as in the (quite different) 
theory mentioned in section [8.7.41 

It turns out (section 15. ip that realizations of our models always have 
jagged routes. A consequence is that (as in Figure 1) a route Tl{^, between 
two points of 5(A) does not pass through any third point ^" of 5(A). This 
prompts the precise definition of geodesic below. 

The route-compatability property is a property that would hold if routes 
were defined as minimum-cost paths, for some reasonable notion of "cost". 
Note that our formal setup does not require routes to be minimum-cost in 
any explicit sense. 

Aspect 2. Subnetworks on locally finite configurations. Here are 
some properties of a fixed feasible subnetwork. 

Lemma 1 Let s be a feasible subnetwork on a locally finite, infinite config- 
uration (Zi). 

(i) The set {r(zj, Zj) n disc(2;, r)}jj- of sub-routes appearing as intersections 
of some route with a fixed disc disc(2;, r) contains only finitely many distinct 
(non-identical) sub-routes. 

(a) For each i and each sequence (zj) with \zj\ oo there is a subsequence 
z'^ = 2;j(fc) and a semi-infinite path vr from Zi in s such that, for each r > 0, 



r{zi,z'^?) ndisc(zi,r) 



vr n disc(2:i, r) for all large k. 
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Proof, (ii) follows from (i) by a compactness argument. To outline (i), if 
false then (by the route-compatability property) the subroutes must meet 
the disc boundary at an infinite number of distinct points, and then (again 
by the route-compatability property) their extensions must meet the bound- 
ary of a slightly larger disc at an infinite number of distinct points, implying 
infinite length and contradicting the "finite length in bounded regions" prop- 
erty dS]) of s. ■ 

Note that Lemma[T]is implicitly about compactness in a topology on the 
space of paths within a given subnetwork s. This is quite different from the 
topology of the space of all subnetworks, mentioned later. 

Terminology: paths, routes and geodesies. A path in s has its usual 
network meaning. Typically there will be many paths between Zi and Zj, 
but (as part of the structure of a feasible subnetwork) one is distinguished as 
the route r{zi, zj). So a route is a path; and a path may or may not be part 
of one or more routes. A singly infinite geodesic in s from Zi is an infinite 
path, starting from Zi, such that any finite portion of the path is a subroute 
of the route r{zi,Zk) for some z^- So Lemma [T]||ii) says that there always 
exists at least one singly infinite geodesic from Zj. A typical point e along a 
route r{zi, Zj) will sometimes be called a path element to distinguish it from 
the endpoints. 

Now write & for the set of all feasible subnetworks s on all locally finite 
configurations x = (xj). It is natural to want to regard 5(A) as a random 
element of ©, which requires specifying a cr- field on &, and as traditional 
we can do this by specifying a complete separable metric space structure on 
& and using the Borel ir-field. 

We outline a "natural" topology in an appendix. In this paper the 
topology plays no explicit role, but one can imagine developments where it 
does - one can imagine constructions using weak convergence, for instance, 
and compactness issues would be key to a proof of the existence part of Open 
Problem [30l However, it might be better to develop such theory within a 
framework where routes are allowed to be rectifiable curves. 

Aspect 3. Desired distributional properties of subnetworks. The 

scale-invariance property (I12|] 

applying Uc to 5(A) gives a network distributed as 5(c~^A) 

is what gives SIRSNs a mathematically interesting structure, and almost all 
our general results in sections [5] and [6] rely on scale-invariance. To indicate 
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how it is used, define £{X) analogously to (fT6]) : 

i{X) := intensity (5 (A)) (21) 

so = I. Then there is a scaling relation 

= A^/2£, < A < oo. (22) 

To derive this relation, consider the scaling map o';^-i/2 that takes 5(1) to 
5(A), and by considering the pre-image A = [0, aV2]2 of the unit square we 
see 

= X area(^) x I 

where the A~^/^ term is length rescaling. 

Similar relations, provable in the same way, will be stated later (j28|301 
[37|) without repeating the proof. 

Aspect 4. Final definition of a SIRSN. Starting from FDDs, a con- 
ceptual and technical issue is how to continue to understand a SIRSN as 
a process over the whole continuum. As an analogy, for continuous-time 
stochastic processes one typically seeks some sample path regularity prop- 
erty such as cadlag. So one might seek some notion of "regularity" for a 
realization, for instance a.e. continuity of routes TZ{zi, Z2) as zi and Z2 vary. 
A version of such continuity is proved, under extra assumptions, in section 
17.21 But as we next explain, in the present context the assumption p{\) < 00 
serves as an alternative regularity condition that enables us to study global 
properties of a SIRSN. 

There are several possible real-world measures of "size" of a road seg- 
ment, quantifying the minor road to major road spectrum - e.g. number of 
lanes; level in a highway classification system; traffic volume. What about 
within our model of a SIRSN? Recalling the definition ()17p of £{X,r), the 
limit 

£{00, r) := UA<oo'f(A,r) 

has (because Ua<ooS(A) is dense) the interpretation of "the set of path 
elements e that are on some route TZ{zi,Z2) with both zi and Z2 at distance 
> r from e" . As shown in section [U assumption (I20p implies that the edge- 
intensity p{r) of £{00, r) is finite and scales as p{r) = p{l)/r. Moreover the 
random process £{00, r) is independent of the sampling process (H(A),0 < 
A < 00) and is an intrinsic part of the global structure of the SIRSN. So 
if we intuitively interpret <£'(oo,r) as "the roads of size > r", then we have 
a mathematically convenient notion of "size of a road segment" emerging 
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from our setup without explicit design. Intuitively, one could view the limit 
£{oo,0+) := Ur>o £{oo,r) as the continuum network of interest. But at 
a technical level it is not clear what are the properties of a realization of 
iS(oo,0+), and we do not study it in this paper. 

3 The binary hierarchy model 

The construction of this model, our basic example of a SIRSN, occupies all 
of section O in several steps. 

• A construction on the integer lattice (sections 13.11 - 13. 4p 

• Extension to the plane (sections 13.51 - 13. 6p 

• Further randomization to obtain invariance properties (section 13. 7p . 
3.1 Routes on the lattice 

For an integer x ^ 0, write height(x) for the largest j E Z+ such that 2-' 
divides x; in other words the unique j such that x = {2k + 1)2^ for some 
/c S Z. Set height (0) = oo. For later use note that in one dimension, any 
integer interval [mi, 7712] contains a unique integer of maximal height, which 
we call peak[mi,m2]. For instance peak[67, 99] = 96 and peak[34, 59] = 48. 

Until section 13.51 we will work on the integer lattice Z^, with vertices 
z = (x, y) whose coordinates have heights > 0. While we are working on the 
lattice it is convenient to use distance ||-Z2 — -Zi||i := |a;2 — + 1^2 — 

Note also that until section \3M we work with deterministic constructions. 

(X) (Y) 

Write Lx and Ly for the lines through {{x,y),y G Z} and {{x,y),x € 

(X) 

Z}. The height of a line Lx is the height of x. 

Fix a parameter 1/2 < 7 < 1. Associate with lines at height h a cost- 
per-unit length equal to j^. Now each path in the lattice has a cost, being 
the sum of the edge costs. Visualize a road network in which one can travel 
along a height-/i road at speed I/7''; so the cost equals time taken. 

Define the route r{zi,Z2) to be a minimum-cost path between zi and Z2- 
There is a uniqueness issue: for instance, for any minimum-cost path from 
(i, i) to {j, j) there is an equal cost path obtained by reflection (x, y) — {y, x). 
However, the estimates from here through section [3^31 hold when r(zi,Z2) is 
any choice of minimum-cost path. We will deal with uniqueness in section 
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A key point of the construction is that if we scale space by 2 then the 
scaled structure on the even lattice (2Z)^ agrees with the original substruc- 
ture on the even lattice, up to a constant multiplicative factor in edge-costs, 
and so the route between two even points will be the same whether we work 
in 1? or (2Z)^. So this "invariance under scaling by 2" property is built into 
the model at the start. 

The fact that moving along the axes has zero cost may seem worrrying 
but actually causes no difficulty (we will later apply a random translation, 
and the original axes do not appear in the final process). Note that the 
cost associated with the line segment from (2'^, 2^*) to (2^*, 0) is 7'*2'* and the 
constraint 7 > 1/2 is needed to make this cost increase with h. Intuitively, 
if 7 is near 1 then the route r{zi,Z2) will stay inside or near the rectangle 
with opposite corners zi,Z2, whereas if 7 is near 1/2 then the route may go 
far away from the rectangle to exploit high-speed roads. 

For this model we will show a property stronger than (|15p : the ratio of 
route-length to distance is uniformly bounded. 

Proposition 2 There is a constant < 00 such that 



Some intuition about possible paths in this model is provided by Figure 
3 (the reader should imagine the ratios of longer/shorter edge lengths as 
larger than drawn). We might have a route as shown in the figure, where 
the two long edges are very fast freeways. But such a route is not possible 
if the fast freeways are too far from the start and destination points. The 
latter assertion will follow from Lemma HI 



Figure 3. Routes like this are possible. 

One might expect some explicit algorithmic description of routes r (21,2:2) 
that one can use to prove the results in sections 13.21 - 13.61 but we have 



len 'c{zi,Z2) < K^\\z2 - zi\\i, ^Z\,Z2 G '2? ■ 
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been unable to do so. Instead our proofs rely on finding internal structural 
properties that routes must have. 

3.2 Analysis of routes in the deterministic model 

Consider the route from zi = {xi,yi) to Z2 = {x2,y2)- The x-values taken 
on the route form some interval Ix 2 [min(xi, 0:2), max(a;i, 0:2)], and simi- 
larly the y- values form some interval ly. Consider the point z* = {x*,y*) = 
{pedk{Ix),pedk{Iy)) and call this point peak*^^V(2;i, 2:2). The notation re- 
minds us that peak^^V(zi, ^2) depends on the route r{zi,Z2), which may not 
be unique. 

Lemma 3 Consider the route r{zi, Z2) from zi to Z2- 
(i) The route passes through z* = peak*^^'*r(2;i, 2:2). 

(a) The route meets the line L^, in either the single point z* or in one line 
segment containing z* (and similarly for ). 

(Hi) Suppose the route passes through a point {x*, y) (for some y 7^ y*) and 
through a point {x,y*) (for some x 7^ x*). Then the route between those 
points is the two-segment route via z* . 

(iv) Suppose z* = z\. If Z2 is in a certain quadrant relative to zi, for 
instance the quadrant [a;i,oo) x [yi,oo), then the route from Zi to Z2 stays 
in that quadrant. 

Proof. We first prove (iii). It is enough to prove that, amongst routes 

between {x*,y) and {x,y*), the two-segment route via z* is the unique 

(Y) 

minimum-cost route. In order to get from {x* ,y) to the line Ly* the route 
must use at least \y — y*\ vertical unit edges; by definition of x* = peak(/3;), 
if these edges are not precisely the line segment from (x*, y) to z* then the 
cost of these edges will be strictly larger; and similarly for horizontal edges. 
This establishes the uniqueness assertion above, and hence (iii). 

For (i), if the hypothesis of (iii) fails then the route must go through 
z*, whereas if it holds then the conclusion of (iii) implies the route goes 
through z* . For (ii), if false then the the route passes through some two 
points {x*,y') and {x*,y") but not the intervening points on that line. But 
(as in the argument for (iii)) the minimum cost path between those two 
points is the direct line between them. 

Finally, (iv) follows from (ii), because if (iv) fails then the route meets 
one boundary of the quadrant in more than one segment. ■ 

Lemma 4 The route r{zi,Z2) from zi to Z2 stays within the square of side 
K'-y 11-22 — -21 111 centered at z\, where K' depends only on 7. 
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Proof. Choose the integer h such that 

2^-^ < \\z2 - 2i||i < 2^. 

As illustrated in Figure 4, there is a square of the form S = [{i - 1)2'^, (i + 
1)2^] X [{j — 1)2^, (j + 1)2''] containing both zi and Z2 (note here i and 
j may be even or odd). We may suppose the route does not stay within 
5 (otherwise the result is trivial). For any point z outside S, call the L°° 
distance from z to S, that is the number d such that z is on the boundary of 
the concentric square Sd = [{i-l)2''-d, {i+l)2^+d] x [{j-l)2''-d, {j+l)2^+ 
d], the displacement of z. Now let d be the maximum displacement along 
the route r{zi, Z2), and choose a point z' along the route with displacement 
d. So the route stays within Sd- 



z' 
• — 



d 


1 


Zi» 






• 

Z2 



Sd 



Figure 4. Construction for proof of Lemma S) 

We may assume, as in Figure 4, that z' is on the top edge of Sd- The 
route needs to cover the vertical distance d between the top edges of S 
and Sd twice (up and down) while staying within Sd, which has side-length 
2/1+1 _|_ 2^_ Now within any integer interval of length a the second-largest 
height H satisfies 2^ < a. So the cost (C, say) of the route outside S is at 
least the cost associated with this second-largest height, which is given by 

dj" where 2^ < 2^+^ + 2d- 

Setting d = b2^, this inequality implies 

log2 C > log2 b + h + {h + l + log2(l + b)) log2 7- 

But for this to be the minimum-cost path, the cost outside S must be less 
than the cost of going round the boundary of S, which is at most x 2^^^"^. 
So 

log2C < hlog2 7 + h + 2. 
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This inequalities combine to show 



log2 6 + (1 + log2(l + b)) log2 7 < 2 

which, because 7 > 1/2, imphes that b is bounded by some constant b^. m 

Lemma U] makes Proposition [2] look very plausible, but to prove it we 
need to extend Lemma[3]to develop internal structural properties that routes 
must have. 

Call a sequence of integers ii, 12, • • • , im a height-monotone sequence from 
ii to im if 

(i) height(ii) > height(22) > . . . > height(im) > 0; 
(h) - ij\ < 2^<='s^^*fe), l<j<m. 

Suppose, for integers mi,m2,m* , we are given a height-monotone sequence 
m* = ii,i2, . . . ,im = "i2 and a height-monotone sequence m* = ji, 32-, ■ ■ ■ ,jq = 
nil. Then we can form the concatenation 

"ii = jqijq-ii ■ ■ ■ )J2,"^*) ^2, ■ ■ ■ ,im = ^2- Call a Sequence that arises this 
way an admissable sequence from mi to m2- See Figure 5. 

5 96 
4 80 
3 72 

2 100 
1 74 

75 99 

(height) 

Figure 5. An admissable path from 75 to 99. This path has range 100 — 72 = 

28. 

Regard a height-monotone or admissable sequence as a path of steps where 
a step from i to j has length \ j — i\. It is clear from (ii) that the length of 
the path in (i) is at most twice the length of the first step. We deduce the 
following crude bound. 

(*) The total length of an admissable path is at most 4 times 
the range of the path, where the range is the difference between 
the maximum and minimum integer points visited by the path. 

Proposition [2] follows immediately from Lemma HI the bound (*) above and 
the following lemma. 
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Lemma 5 The route from z\ = {xi,yi) to Z2 = {x2,y2) consists of alternat- 
ing horizontal and vertical segments, in which the successive distinct x-values 
of the segment ends (the turning points) form an admissable sequence from 
xi to X2, and the successive distinct y-values form an admissable sequence 
from yi to y2- 

Proof. In view of Lemma[3]we can reduce to the case where zi = peak^^V(2:i, Z2), 
and we need to show that the successive distinct x-values form a height- 
monotone sequence, as do the y-values. Without loss of generality suppose 
that xi < X2, that yi < y2 and that the first segment is horizontal. So the 
route is of the form 

(.xi,yi) = (a;(i),y(i)) (x(2),y(i)) ^ (.X(^2),y{2)) ^ {x{3),y(2)) ^ • • • 

It suffices to show that for each edge of the route, say the edge — )■ 
and for each point (say on the edge other than the 

starting point, we have height(x*) < height(x(j)). This is true for the first 

two edges of the route by definition of zi as peak''^V(zi, 22). Suppose it fails 
first at some point (x*,y(j)). Then the route has proceeded y(i-i)) 
(x(j),y(j)) — (x*,y(j)) instead of the alternate path via (x*, Now 
inductively height < height so the cost of the horizontal edge 
is less in the alternate path; so for the route to have smaller cost it must 
happen that the cost of its vertical edge is smaller than in the alternate 
path, that is height(x(j)) > height(x*), contradicting the supposed failure. 
■ 

3.3 Further technical estimates 

The next lemma will be key to bounding network length, more specifically 
to showing ^ < 00 later. 

Lemma 6 There exists an integer b > 1, depending only on 7, such that for 
allh>0 and all rectangles of the form [^2''+^ {i+l)2^+'']x[j2^ , {j+l)2^], the 
route r(zi,Z2) between two points zi,Z2 G outside (or on the boundary of) 
the rectangle does not use any horizontal edge strictly inside the rectangle. 

Note there may be routes using a vertical line straight through the rectangle. 
Proof. Suppose false; then there are two points zi,Z2 on the boundary 
of the rectangle such that the route between them lies strictly within the 
rectangle and contains a horizontal edge. Because the speed on an interior 
edge is less than the speed on a parallel boundary edge, this cannot happen 
when zi and Z2 are in the same or adjacent boundaries of the rectangle. 
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because the path around the boundary is faster. Suppose they are on the 
top and the bottom boundaries. Then the height of the horizontal edge is less 
than the heights of the starting and ending y-values, contradicting Lemma 
[5l The only remaining case is when zi and Z2 are on the left and right 
boundaries. Using Lemma [S] again, the route cannot use a vertical edge 
inside the rectangle, so the only possibility is a single horizontal segment 
passing through the rectangle. Such a path has cost at least 2'*"'"* t'*"^) 
because the height of the line is at most h — 1, whereas the path around 
the boundary has cost at most 2'^"'"* 7'^ + 2^ 7^+''. So the potential route is 
impossible when 2^ + 7* < 2*7"^ which holds for sufficiently large b. m 

Corollary 7 If a route r(zi,Z2) uses a height-h segment through zq, then 
mindl^i — zqIIi) 1 1 -22 — -^oHi) ^ 2''(2^ + 1) for b as in Lemma\^ 

Proof. Consider a unit-length horizontal (without loss of generality) edge of 
height h at zq. It is in the interior of some rectangle of the form [z2'*"'"^"'"*, {i + 
1)2^+1+''] X [i2'*+i, (j + l)2'*+i]. By Lemma [6] applied with /i + 1, either zi 
or Z2 must be within that rectangle. ■ 

Perhaps surprisingly, we do not make much explicit use of the determin- 
istic function cost (zi, 2:2) giving the cost of the minimum-cost route in the 
integer lattice, but will need the following bound. 

Lemma 8 There exists a constant K" such that 

cost(zi, Z2) < K"\\z2 - Zl\f 
where (3 := log(27)/log 2. 

Proof. As in Figure 4 in the proof of LemmaSl there is a square of the form 
S = [{i- 1)2^, {i + 1)2^^] X [(j - 1)2'', {j + 1)2^] containing both zi and Z2, 
where h is the integer such that 2^^^ < ||z2 — ziHi < 2'^. As observed there, 
the cost of going all around the boundary of S is 0{j^2^). By considering 
a path from zi using the "greedy" rule of always switching to an orthogonal 
line of greater height, it is easy to check that the cost of this greedy path 
from zi to the boundary of S is also 0{'j^2^). Hence cost(zi, ^2) = 0{'y^2^) 
and the result follows. 

3.4 Finessing uniqueness by secondary randomization 

As previously observed, minimum-cost paths are not always unique. We 
conjecture that, at least when 7 is not algebraic, there is some simple classi- 
fication of when and how non-uniqueness occurs. But instead of addressing 
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that issue we can finesse it by introducing randomness (which we need later, 
anyway) at this stage. One possible way to do so would be to use the uniform 
distribution on minimum-cost paths. Instead we use what we wih cah sec- 
ondary randomization to choose between non-unique minimum-cost paths. 
Place i.i.d. Normal(0, 1) random variables ("weights") Ce on the edges e of 
7?. Any path has a weight Yle in path "^e- Define the route Tlo{zi, Z2) to be 
the minimum- weight path in the set of minimum-cost paths from z\ to 2:2- 

3.5 Extension to the binary rational lattice 

The notion of height extends to binary rationals: if x G M is a binary rational 
and x 7^ 0, write height (x) for the largest j £ I' such that 2^ divides x; in 
other words the unique j such that x = {2k + 1)2^ for some /c G Z. 

For —00 < < 00 let be the lattice on vertex-set {2^ z : z G Z^}, 
in other words on the set of points in whose coordinates have height 
> H. So far we have been working on the integer lattice Z^, but now the 
results we have proved extend by (binary) scaling to analogous results on 
the lattices Z|^. We will use such scaled results as needed. 

Note in particular the following consistency condition as H varies. Take 
Hi < H2. Consider the route, in Z|^^, between two vertices of Z^^. By 
Lemma [5] and the definition of admissable, any minimum-cost path stays 
within the lattice Z|^^ . So the set of minimum-cost paths is the same whether 
we work in Z|^^ or in Z^^. Note also that each edge e in Z|^ corresponds 
to two edges 61,62 of Z^_-|^. So we can couple the edge- weights by making 
Ce = Cei +Ce2 (ouly this infinite divisibility property of the Normal is relevant 
to the construction) and this gives a "consistency of secondary weights" 
property, which implies that the random route TZq{zi,Z2) is also the same 
whether we work in Z|^^ or in Z^^. 

So we have now defined random routes Tlo{zi, Z2) for all unordered pairs 
of vertices in Z^^^ := U/f>_ooZ|^. From the "minimality" in the construc- 
tion it is clear that the routes satisfy the route-compatability properties 
(iii,iv) from section 12.21 

3.6 Extension to the plane 

We want to define routes TZq{zi,Z2) between general points zi,Z2 of as 
H — —00 limits of the routes TZq{z^ ^ Z2) between vertices such that 

GZ^,, zf ^z, (i = l,2) (23) 
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Proposition [9] formalizes this idea. The proof in this section is the most 
intricate part of the construction, which can thereafter be completed (section 
13. 7p by "soft" arguments. 

As a first issue, what does it mean to say that, under (|23p . 

routes r{z^,Z2 ) converge to a route r{zi,Z2)'^ (24) 

We define this to mean: 

for each Hq > — oo, the subroute yhq{zi , Z2) consisting of path 
segments of r(z^,2;|^) within lines of height > Hq is, for suffi- 
ciently large negative H, a path not depending on H - call this 
path rHo(zi,Z2). 

When this property holds. Lemma [5] implies that rHo{zi,Z2) is a connected 
path, consistent as Hq decreases, and using Proposition [2] and scaling we see 
that the closure of Uho>-oo'''Ho{zi, Z2) defines a route r{zi, Z2) satisfying the 
"jagged" condition of section \T2[ 

Proposition 9 There exists a subset A C of zero area such that, if zi 
and Z2 are outside A, there exists a random jagged route 1Zq{zi^Z2) such 
that, whenever [2^) holds, then (24\ ) holds. 

The proof relies on the fact that, for particular configurations illustrated in 
Figure 6, routes from a certain neighborhood to distant destinations must 
pass through a particular point. In fact all that matters is the existence 
of such a configuration, not the particular one we now exhibit. Consider a 
square G= [2'^,2'^ + 2]2 and points 6 = (2'^ + l,2'^) and d = (2'^ + l, 2'^+2-''), 
illustrated in Figure 6. What is relevant is the heights of the lines involved, 
indicated in the figure. 



22 



NW 



b* 



c* 



1 



NE 







w 



E 



-h 



h 



sw 



b 



SE 



h 







1 



Figure 6. The big square G and the small square S. Marginal labels attached 
to lines are line-heights, not coordinates. 

Lemma 10 There exist large h and small e (depending on such that, 

in the configuration shown in Figure 6, every route from inside the small 
square E := d + [— £, 0] x [0, e] to the boundary of G passes via b. 

Proof. For each point c on the boundary of G there is a counter-clockwise 
path 7ri(6, c) and a clockwise path TT2{b,c) along the boundary from b to c. 
These paths have equal cost for the point c* = {2^ + 2 — 7'^"^ , 2^* + 2), which 
is near the NE corner point of G. We will need the following lemma. 

Lemma 11 There exists h such that the following hold. 

(a) The paths 7ri(6, c*) and 712(6,0*) attain the minimum cost over all paths 
from b to c* , and are the only paths to do so. 

(b) The only minimum-cost paths from d to c* are the two paths consisting 
of the segment [d, b\ and the paths 7ri(5, c*) or 7r2(6, c*). 

(c) There exists 77 > such that any path from d to c* that avoids the 
segment [d, b\ has cost at least rj greater that the minimum- cost paths. 

Note a technical point. We are working on := Uh>-oo^h ^^'^ ^ 
not be in Z?.^. To be precise we should replace c* in the arguments below 
by a sequence — >■ c*, but that requires awkward notation we prefer to 
avoid. 
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Granted Lemma [TT] we deduce Lemma [TOl as follows. Consider a point 
c on the counter-clockwise path from b to c* (the clockwise case is similar). 
Then the following must hold, because any counter-example path to c could 
be extended along the boundary from c to c* and would give a counter- 
example to Lemma [TTl 

(a) The path 7ri(6, c) is the unique minimum-cost path from b to c. 

(b) The path consisting of the segment [d,b] and the path 7ri(5, c) is the 
unique minimum-cost path from d to c. 

(c) Any path from d to c that avoids the segment [d, b] has cost at least r] 
greater that the minimum-cost path. 

Lemma [8] extends by scaling to Uh>-oo'^'h, and so the function cost(-,-) 
extends to a continuous function on M^. So we can choose H so that the 
square S = d+[-2-^,0] X [0,2-^] satisfies 

It is easy to check that a minimum-cost path from s £ T, to d does not 
meet [d, b] except at d. 

Consider s E E and a point c as above. So cost(s, c) < cost((i, c) + rj/3. 
Suppose a minimum-cost path from s to c does not meet the segment [d, b] . 
Then the path from d to c via s would have cost < cost(d, c) + 2r//3 and 
would not meet [d,b], contradicting (c). So a minimum-cost path from s to 
c must meet the segment [d,b]. Then, by uniqueness in (b) (for the path 
from d to c), it must continue via b, establishing Lemma [TOl ■ 

Proof of Lemma llll (I thank Justin Salez for completing the details of 
this proof.) In outline, we use the "structure of paths" results in Lemmas [3] 
and [5] to reduce to comparing costs of a finite number of possible routes. We 
will make use of the following preliminary observations, which are straight- 
forward to check : 

(i) the only minimum-cost path from SE to NW is SE — t- SW — t- NW ; 

(ii) the only minimum-cost paths from SW to NE are SW — )• NW — )• NE 
and SW ^ SE ^ NE ; 

(iii) O — 7- 6* — 7- NE is a minimum-cost path from O to NE. 

Consider assertion (a). The cost associated with paths vri(6, c*) and '7T2{b, c*) 
equals 27^ + 27, which (by choosing h large) is less than 2. Now consider 
some minimum-cost path vr from 6 to c*. Since both end-points have their 
y— coordinate at height > 1, all horizontal segments of vr must have height 
> 1 (Lemma [5]). In other words, the length of every vertical segment must 
be an even integer. If the last vertical segment of vr were ending strictly 
between NW and NE, then its cost would be at least 2, contradicting 
optimality. Thus, tt must pass through NW or NE, and observationss (i) 
or (ii) complete the proof. 
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Now consider assertions (b) and (c). Let vr be any path from d to c*, and let 
z be the point at which vr first meets the boundary of the rectangle formed 
by {E^W^SW^SE\. Let 7r',7r" denote the subpaths of vr from d to z and 
from z to c*, respectively. There are four possible cases : 

• z S (i?, VF) : since all segments in vr' have height < 0, replacing vr' by 
d — >• O — 7- z cannot increase the overall cost. In the resulting path, 
one may further replace the subpath from O to c* by O — t- 6* — ?■ c* 
without increasing the cost, by (iii). This shows : 

cost(7r) > 2 + 7 - 2"'' - 7^". 

• z € (VF, : all horizontal segments in vr' have height < —1, so 
cost(7r') > 7~^. By (ii), one also has cost(7r") > cost(2; — > NW — > c*). 
Combining these two facts yields 

cost(7r) > 27 + 7""*". 

• z E {SE, E) : replacing the subpath vr" by z — )• NE — )• c* cannot 
increase the overall cost, by part (a). In the resulting path, the subpath 
from d to E costs at least 7^"*^ + 7(1 — 2"^^), because the horizontal 
and vertical heights are < — 1 and < 1, respectively. Thus, 

cost(7r) > 27 + 7""*". 

• z G {SW,SE) : all segments of vr' have height < — 1 except those 
included in [0,b], which have height 0. Thus, 

cost(7r') - cost((i 6 z) > (7"^ - l)len {[b,d] \ vr') . 

Moreover, by part (a), cost(6 — )■ z) + cost(vr") > cost(vri). Thus, 

cost(vr) > (7^"'^ — l)len {[b,d] \ vr) + cost((i — )• 6) + cost(vri). 

Let us sum up: in the first three cases, the cost of vr exceeds that of our 
two candidates by at least 1, for h sufficiently large. In the fourth case, the 
excess is at least (7"^ — l)len ([6, d] \ vr). This proves both (b) and (c), with 
r] = 2-^{-/-'^ -1). m 

Proof of Proposition [9]. In each basic 2^^^ x 2^^^ square G of Z^^^i there 
is a copy of the Figure 6 configuration; let Tiq be the corresponding small 
square. Let B := Ug'^g be the union of those squares, for the fixed h given 
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by Lemma [TOl Then for i > 1 let Bi := a2-iB be rescalings of B. Each Bi 
has the same density, which by Lemma [10] is non-zero, and a straightforward 
use of the second Borel-Cantehi lemma (with sufficiently well-spaced values 
of i) shows that the set 

j4 := {z G : z in only finitely many B^} 

has area zero. 

Now consider zi G A'^. Then there exists a sequence ij = ij{zi) — t- 
oo such that zi £ Bi^ and the associated hi.{zi) — )• zi. Consider Z2 7^ 
zi and {zi ,Z2) — )■ {zi,Z2) as in (f23|) . For j larger than some jQ{zi,Z2), 
Lemma 1101 implies that for all sufficiently large H the route TZo{z^ , Z2) 
passes through 6,^. (zi ) . But the routes between the b^. {zi),j > 1 are specified 
by the construction on Uh>-oo'^'h- It follows that, when zi and Z2 are both 
in A'^ we have convergence in the sense of (fM|) to a route TZo{zi,Z2)- m 

We digress to give the technical estimate that will show £ < 00 in this 
model. 

Lemma 12 For the routes TZq in Proposition\^ take the union over points 
of a rate-1 Poisson point process S(l) of the routes TZoi^i^'); o,nd let 
S* he the intersection of that union with the interior of a unit square U = 
[i,i -|- 1] X + 1]. Then the expected length of S* is at most 2^^"^ , for b as 
in Lemma\^ 

Proof. Lemma [6] was stated for h>Q and vertices in Z^, but by scaling 
it holds for /i < and vertices in M^. Consider /i < 0. Within U there are 
2-/1-1 horizontal unit-length line segments at height /i, and these can be 
split into 2^^^~^ segments of length 2^. Consider such a line segment, C, say. 
It is in the interior of some rectangle of the form [i2^^^^^ , {i + 1)2''+^+''] x 
[j2^^^ , [j -\- 1)2'*^"^]. By Lemma El applied with h + 1, the only possible way 
that the segment C, can be in a route 'R-oii^i') is if ^ or S,' is within the 
rectangle. (And the same holds for any piece of C,, by considering a sub- 
rectangle). The chance the Poisson process contains such a point is at most 
the area of the rectangle, which is 2^'^+^"^^. 

So the contribution to mean length from a particular segment C, is at 
most 2^ X 2^'^+^+'', and then the contribution from height-/i horizontal lines 
is at most 2^ x 2^^+"^+^ x 2~2/i-i ^ 2^+^+^. Summing over h < -I and 
adding the same contribution from vertical lines gives the bound 2-^"*"''. ■ 
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3.7 Completing the construction by forcing invariance 

Proposition E] gives paths 1Zq{zi,Z2) when zi^Z2 G A^. The process IZq can- 
not be translation- or rotation- invariant (in distribution), because the axes 
play a special role (infinite speed); though by construction the process is in- 
variant under a2 (scaling space by a factor 2). But there is a standard way of 
trying to make translation-invariant random processes out of deterministic 
processes, by taking weak limits of random translations of the original pro- 
cess. In our setting this can be done fairly explicitly as follows. For u G 
let Tu be the translation map Tu{z) = u + z, z G on points, and let 
act on routes in the natural way. Take Un uniform on the square [0,2"]^, 
and couple the random variables {Un,n > 1) by setting Un = Un+i mod 2" 
coordinatewise. Define 

7^(")(Z1,Z2) = T^uMizi + Un,Z2 + Un)). (25) 

In words, translate points by Un, use TZq to define a route between the 
translated points, and then translate back to obtain a route between the 
original points. 

Now the only way that TZ^"'^^\zi, Z2) could be different from Tl^"'\zi, Z2) 
is if the route TZo{zi + Un,Z2 + Un) intersects the boundary of the square 
[0,2"]2, which, using Lemma m has chance 0(2 "). So we can define a 
random network TZi_[ via the a.s. limits 

7^t-i (-21 1 Z2) = Tl^"'\zi, Z2) for all sufficiently large n. (26) 

This process is translation-invariant, because for fixed z G'M? the variation 
distance between the distributions of C/„ and Un + z mod 2" tends to zero. 

For < c < 00 write cfc for the scaling map z — t- cz on M^, and recall 
that TZq is invariant under (T2. Now for TZ^"'^ at (I25p . 

^727^(")(zl, Z2) = a2T_uMzi + Un, Z2 + Un) 

= T_2Un(^2'R'*{zi + Un, Z2 + Un) 

= T_2(7„'^*(2zi + 2Un,2z2 + 2Un) by invariance of TZq under a2 

= T^u„+iT^*{2zi + Un+i,2z2 + Un+i) because Un+i = 2Un 

= 7^("+l)(2zl,2z2). 

Hence the distribution of the limit 7^t-i is invariant under £72- 

Of course our construction so far is not rotationally invariant, but ap- 
plying a uniform random rotation to 7^t-i gives a network TZj-_[ whose dis- 
tribution is invariant under rotation, as well as preserving distributional 
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invariance under translation and under a2 ■ Finally, we get a process TZ with 
scale-invariant distribution by random rescaling via the scale-free distribu- 
tion: 

TZ{zi,Z2) = <yi/cT^r-i {Czi,Cz2), P(C e dc) = ^r^, 1< c < 2. (27) 

This completes the construction of the binary hierarchy model IZ. To 
check it satisfies the formal setup of a SIRSN in section 12.21 the only re- 
maining issue is to check that the parameters EDi,^ and p{l) are finite. 
For the former, Proposition [2] implies the corresponding bound in terms of 
Euclidean distance 

len TZq{zi,Z2) < 2^/'^K^\\z2 - zi\\2 

and this bound is unaffected by the transformations taking IZq to IZ. So 
EDi < 2^/'^K^. For I, in the notation of Lemma [T2l the edge- intensity 
of Ug^5/gH(i)'^o(C,'^') is at most 2''+2 + 2, the "+2" terms arising from the 
edges of 1?. This edge- intensity is unaffected by the transformations taking 
7^0 to 7^i-i • Scaling by C in (j27|) multiplies edge-intensity by C, so finally 
I < (2''+2 + 2)EC. To bound set r{h) = 2^{2^ + l). Corollary [7] implies 
that, for routes IZq, if an edge element is in a route between some two points 
at distance > r[h) from the element, then the edge has height > h. The 
edge-intensity of edges with height > h equals 2^~^. These quantities are 
unaffected by translation and rotation; and the scaling by ctc can at most 
increase the edge-intensity by 4. So the edge-intensity in IZ of £{\,r{h)) is 
p{X,r{h)) < 4 • 2^~^ Choosing h such that r{h) < 1 we deduce p{l) < oo. 

3.8 Remarks on section [31 

The "combinatorial" arguments in sections 13.11 - 13.41 are obviously specific 
to this model. But the property implicit in Lemma fTOl (that there exist con- 
figurations in which all long routes from a small neighborhood exit the unit 
disc at the same point) is closely related to desirable structural properties 
of SIRSNs discussed in section [3 

Lemma [20] later shows that in general i < 2p{l), so our argument above 
that ^ < oo could be omitted, though it is pleasant to have a self-contained 
construction. 

4 Other possible constructions 

The model in section 3 has some very special features, in particular that in 
any realization we see a (scaled and rotated) square lattice of roads. Below 
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we outline two other constructions which, we conjecture, produce SIRSNs, 
the technical dificulty being to prove a.s. uniqueness of routes defined as 
minimum-cost paths. 

4.1 The Poisson line process model 

For each m = 1,2,3,... take a rate-1 Poisson line process, and attach 
Uniform(m — l,m) marks to the lines; the union of all these is a Poisson 
line process with "mark measure" being Lebesgue measure on (0, oo). By a 
one-to-one mapping of marks one can transform to the mark measure with 
density x~"' on < x < oo, where we take the parameter 2 < 7 < oo. So in 
any small disc, there is some finite largest mark amongst lines intersecting 
the disc. 

Picturing the lines as freeways and the marks as speeds, for any pair of 
points zi, Z2 on the lines there is some finite minimum time t{zi, Z2) over all 
routes from zi to Z2, and analogous to Lemma[8]one can show (Wilf Kendall: 
personal communication) that this function extends to a random continuous 
function i[zi,Z2) on the plane. The technical difficulty is to show that for 
given {zi,Z2) there is an a.s. unique route attaining that time; if that were 
proved, establishing the remaining properties required of an SIRSN would 
be straightforward. In particular, scale-invariance would follow from the 
form x~'^ of the mark density. 

4.2 A dynamic proximity graph model 

This potential construction of a SIRSN is based on a space-time Poisson 
point process (H(A),0 < A < 00). Note that to study such a SIRSN one 
would use an independent Poisson point process to define 5(A). Note also 
that the corresponding "static" model, called the Gabriel network, is a mem- 
ber of the family of proximity graphs described in jlSJfS]; any family member 
could be used in the construction below. 
Here's the construction rule. 

When a point arrives at time A, consider in turn each existing 
point ^' G H(A— ), and create an edge (C,^') if the disc with 
diameter {£,,£,') contains no other point of H(A— ). 

Write G(A) for the time-A network on points H(A). Note the automatic 
scale-invariance property 

the action of cJc on G(A) gives a network distributed as G(c~^A). 
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Now fix a parameter < 7 < 7* for some sufficiently small 7* > and view 
an edge created at time A as a road with speed A^'*'. Defining routes in G(A) 
as minimum-time paths, it seems intuitively plausible, as in the Poisson 
line process model, that that we can extend the minimum-time function on 
UaH(A) to a continuous function t{zi,Z2) and then prove there is an a.s. 
unique route attaining that time. Again, if that were proved, establishing 
the remaining properties required of an SIRSN would be straightforward. In 
particular, scale-invariance would follow from the fact that the construction 
rule is scale-invariant. 

5 Properties of weak SIRSNs 

In this section we study properties that hold for any weak SIRSN, that is 
when we do not require ([20]) but instead require (fT6|) . These are essentially 
properties of the sampled subnetworks 5(A) for fixed A - we cannot get 
A — )• 00 results. 

5.1 No straight edges at typical points 

If a point ^ of H(A) is the start of some straight line segment of length > r 
in 5(A) then consider the subroutes of length exactly r from ^. The edge 
process of such subroutes has some edge-intensity i(A,r). In independent 
copies of H(l) these edge-processes cannot have any positive-length overlap. 
So by regarding H(n) as the union of n independent copies of H(l) we have 
L{n,r) = ni(l,r). But by the general scaling property 

L{X,r) = X^/h{l,rX^/^). (28) 

Since /-(l,r) < ^ < 00 these two different scaling relations imply i(l,r) = 
for all r > 0. 

This proves (a) below; note the consequence (b), implied by the definition 
of feasible path in the section \T2\ setup. 

Proposition 13 S{X) has the following properties a.s. 

(a) 5(A) contains no line segment [S,,z] of positive length, for any ^ G H(A). 

(b) The route '7^(Ci,^2) between two points 0/ H(A) does not pass through 
any third point ^3 o/H(A). 

5.2 Singly and doubly infinite geodesies 

Recall from section [23] that a singly infinite geodesic from a point ^0 in 5(A) 
is an infinite path, starting from such that any finite portion of the path 
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is a subroute of some route 71(^0, Lemma [T] showed 

There is a.s. at least one singly infinite geodesic from each point of 5(A). 

(29) 

A doubly infinite geodesic in 5(A) is a path vr which is an increasing union 
of segments vTfc, where each vrfc is a segment of some route 7^(^fc,^^) between 
two points of H(A), and both endpoints of vr^ go to infinity. 

Previous work on very different (e.g. percolation- type |19j ) networks 
suggests there may be a general principle: 

In natural models of random networks on M? or Z^, doubly infi- 
nite geodesies do not exist. 

Proposition [T3] proves this for weak SIRSNs based on a simple scaling argu- 
ment. Note however this argument depends implicitly upon our assumption 
^ < oo which seems rather special to our setting. 

Recall the setup of ([TTlfTS]) . S{X, r) C 5(A) is the set of points z in edges 
of 5(A) such that z is in the route TZiCjC) for some of H(A) such that 
min(|z — ^1, \z — ^'D > r. And p(A,r) is the edge-intensity of £{X,r). By 
scaling, 

p{X,r) = X^/'^p{l,rX^/^) (30) 

Proposition 14 p{X,r) — )• as r — )• oo. In particular, 5(A) has a.s. no 
doubly infinite geodesies. 

Proof. For fixed A the edge-processes £{X,r) can only decrease as r in- 
creases, and the limit £{X, oo) := rir£{X, r) is by definition the set of path el- 
ements in doubly infinite geodesies. This limit has edge-intensity p{X, oo) = 
limr_).ooP(A, r) > 0. So it is enough to prove |?(A,oo) = 0. Suppose not. 
Then by the scaling relation (I30p 

p{X, oo) = A^/^p(l, oo), < A < oo. 

We claim that in fact 

£{X, oo) = £{l,oo) a.s. for A < 1, 

which (because we know p(l,oo) < oo) implies p(l,oo) = and completes 
the proof. 

To prove the claim, note that for any finite-length segment ttq of a doubly 
infinite geodesic in 5(1), there are an infinite number of distinct pairs CjjCj 
of H(l) such that Tl{Cj,C'j) contains ttq, and for each pair there is chance A^ 
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that both points are in H(A). These events are independent (because H(A) 
is obtained from H(l) by independent samphng) so a.s. an infinite number 
of pairs are in H(A), implying that ttq is in a doubly infinite geodesic 

of 5(A). ■ 

Remark. The limit used here is different from the limit p{r) := lim\^^ p(A, r) 
featuring in assumption (j20p . 

5.3 Marginal interpretation of i 

Recall i is defined as the edge- intensity of 5(1), which is the subnetwork 
on a rate-1 Poisson point process H(l). Now augment the network 5(1) 
by including the point at the origin and the routes from the origin to each 
^ G "(I)- The newly added edges have some random total length L. 

Proposition 15 EL = £/2. 

Proof. Recall ([22]) the scaling relation ^(A) = X^^'^i, where £{X) is the edge- 
intensity of the subnetwork 5(A) of a Poisson process of point-intensity A . 
Differentiating with respect to A, 

/(I) = It 

So we need to show EL = £'{1). 

Consider the space-time Poisson point process (H(A),0 < A < oo) from 
section [2^21 Each arriving point creates some additional network length, say 
L(,^), and for a point arriving at time A, write £{X) for the mean additional 
network length. Now 

i{xo) = E Yi ^(^) = r ^~(^) 

€eS(Ao)n[o,i]2 "^o 

and so f (1) = £{1). 

5.4 A lower bound on network length 

Write A for the parameter EL^i of a SIRSN. Write ^^(A) for the minimum 
possible value of £ in a SIRSN with a given value of A. 

Proposition 16 £*(A) = n{{A - 1)"^/^) as A | 1. 
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The proof is based on a bound (Proposition [T7|) involving the geometry 
of deterministic paths, somewhat similar to bounds used in [8] section 4. 
Figure 7 illustrates the argument to be used. 
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Figure 7. 

Proposition 17 Let a,L,D and Bq be positive reals satisfying 9q < Tr/2 
and 

2^(Z)-l)2 + (fL + l)2 = (l + 2a)(3L + 2) (31) 

^(^-l) = 4«V(3i + + 1. (32) 

-Let TZ be a route from some point z\ in the unit square [—1,0] x [0,1] to 
some point Z2 in the unit square [3L, 3L + 1] x [0, 1], and suppose 

len(7^) < (1 + 2q)|22 (33) 

Take U uniform random on [L,2L]. The route TZ first crosses the vertical 
line {{U,y), —oo < y < oo} at some random point {U,(,{U)) and at some 
angle f3{U) £ (— ^, ^) relative to horizontal. Then 
(^)\m\<D. 

(ii)n\i3m<0o)>i 

Proof. The maximum possible value of S,{U) arises in the case where zi = 
(— 1, 1), Z2 = (3L + 1, 1), U = |L, the route consists of straight lines from 
zi to {U,S,{U)) to Z2, and the route- length attains equality in ([33|) . In this 
case the value of ^(C/) is the quantity D satisfying (f3T]) . establishing (i). 

Writing /3{u) for the angle (relative to horizontal) of the route at x- 
coordinate n, then the length (A, say) of the route between x-coordinates L 
and 2L equals J^^ cos]3(u) This implies 

A-L>i^-l)xL¥{l3{U)>eo). 
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But by considering excess length (relative to a horizontal route) , ()33p implies 

A - L < 2a\z2 -zi\< 2aV(3L + 2)2 + 1. 

Combining these inequalities gives a lower bound on P(/3([/) > ^o) which 
equals 1/2 when satisfies (f32]) . establishing (ii). ■ 

Proof of Proposition 1161 Consider a SIRSN with parameters £ and A 
and with induced subnetwork 5 on a Poisson point process H. Set a = 
A — 1. Suppose we can choose L,D,6q to satisfy, along with the given a, 
the equalities (131|32p ~ note this leaves us one degree of freedom. 

With probability (1 — e"^)^ there are points zi and Z2 of the Poisson 
process in the unit squares [—1,0] x [0,1] and [3L,3L + 1] x [0,1]. By 
Markov's inequality and the definition of A, with probability at least 1/2 
the route TZ{zi, Z2) has length at most (l + 2a)|2;2 — -^il- Applying Proposition 
[T7]we deduce that, with probability > (1 — e~^)^/4, the network S contains 
an edge that crosses the random vertical line {{U,y) : —00 < y < 00} 
at some point {U,^{U)) with —D < S,{U) < D and crosses at some angle 
(3{U) S (— 6*0,6*0) relative to horizontal. 

If we translate vertically by 2D, to consider routes between the unit 
squares [-1, 0] x [2D, 2D+1] and [3L, 3L+1] x [2D, 2D+1], then the potential 
crossing points (using the same r.v. U) for the translated and untranslated 
cases are distinct. Now by considering translates by all multiples of 2D, 
and noting that the distribution of crossings of the random vertical line 
{{U, y) : —00 < y < 00} is the same as for the y-axis, we have shown 

the mean intensity of crossings of the network S over the y-axis 
at angles G {—6q,6()) relative to horizontal is at least ^^"g^^-, ^ . 

The stochastic geometry identities ([21 E]) relates this mean intensity to the 
parameter i via 




Combining with the previous inequality we find 



- 21D6lo' 



Now set L = a and consider the solutions of (|31l32p in the limit as 
a 4, 0: we find that solutions exist with 



6*0 ~ V24^; D^IO 



which establishes Proposition 1161 
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5.5 The minimum value of I and the Steiner tree constant 



Take k uniform random points Zi, . . . , in a square of area k and con- 
sider the length LsT{k) of the Steiner tree (the minimum-length connected 
network) on Zi , . . . , Z^. . Well-known subadditivity arguments [2T| B5j imply 
that KLsrik) ~ Cst^ for some constant < Cst < oo. One can define Cst 
equivalently (see [3] for results of this kind) as the infimum of c such that 
there exists a translation-invariant connected random network over H(l) 
with edge-intensity c. From the latter description it is obvious that in any 
SIRSN we have i > Cst- So the overall infimum 

:= infimum of £ over all SIRSNs (34) 

satisfies > Cst, and below we outline an argument that the inequality is 
strict. First we derive some simple lower bounds on Cst and 

(i) Write 6(^) for the distance from ^ to its closest neighbor in H(l). The 
discs of center ^ and radius h{^)/2 are disjoint as ^ varies and must contain 
network length at least 6(^)/2, so 

Cst > i]E6(0 = \. 

(ii) In a network of edge-intensity c, ([2]) shows the mean number of edges 
crossing circle(0, r) equals 2-Kr x 2-k~^c = Arc. If there is a point of H(l) 
inside disc(0, r) then there must be some such crossing edge, so 

1 — exp(— vrr^) < 4rc. 

So 

1 — exp(— vrr^) 
Cst > sup ^ ^ 0.283. 

r 4r 

(iii) We can get a better bound on i^, by using Proposition [15] as follows. 
Using the intensity calculation above, in a network of edge-intensity i the 
probability that no edge crosses circle (0,r) is at least 1 — 4r£. When a new 
point arives at ^ in the 5(A) process at time A = 1, if no existing edges cross 
circle(^, r) then the added network length L is at least r. So 

EL > supr(l -4r^) = j^^. 

r 

But Proposition [T5] says I = 2EL and so we have shown 

4 > \/l78 ~ 0.353. (35) 
One could no doubt obtain small improvements by similar arguments. 
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Here is an outline argument that > Cqt- 

(i) In the Steiner tree on the Posson point process vertices of degree 
> 1 have non-zero density, and their edges meet at some varying angles, 
whereas at the Steiner points (non-vertex junctions) edges must meet at 
120 degree angles. 

(ii) If there were a SIRSN with i ~ Cst, then 5(1) would have essentially the 
properties (i). But then in 5(1/2), obtained by deleting half the vertices of 
H(l) to get H(l/2), some of the deleted vertices would remain as junction 
points. The "varying angles" property implies the edge-intensity -^(1/2) of 
5(1/2) is strictly larger than that of the Steiner tree on H(l/2), contradicting 
the scale- invariance property that the edge-intensities of 5(1) and the Steiner 
tree on H(l) are essentially equal. 

6 General SIRSNs and their properties 

In this section we study properties of 5(A) in the A — >• c« limit, for a general 
SIRSN. Roughly speaking, this is studying "the whole SIRSN" instead of 
sampled subnetworks, and such results depend on assumption (|20]l . 

Recall again the setup from (fT7|) - (pO|) . So p{X,r) is the edge-intensity 
of £{X, r), which is the process of points z in edges of 5(A) such that z is in 
the route 7?.(^,^') for some in H(A) such that min(|z — ^|, \z — ^'|) > r. 
Recall also from (pO]l the scaling relation p(A,r) = A^/^p(l, rA^/^). Defining 

p{r) := lim p{X,r) < oo (36) 

A— )-oo 

the assumption (pOl) that p{l) < oo and scaling imply 

p(r) = p{l) X r^^, < r < oo. (37) 

6.1 A connectivity bound 

Assumption (|20p has a direct implication for the qualitative structure of 
a SIRSN: all the routes linking two regions, once they get away from a 
neighborhood of the regions, use only a finite number of different paths. We 
first give a version of this result in terms of discs. 
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Figure 8. Schematic for routes from inside disc(0, 1/2) to outside disc(0, 3/2) 
crossing the unit circle. 

Proposition 18 Take < r < 1 and let N{X, r) be the number of distinct 
points on the unit circle at which some route R{£,, £,') between some ^ G 
S(A) n disc(0, 1 - r) and some ^' G H(A) n (M^ \ disc(0, 1 + r)) crosses the 
unit circle. Then 

E lim N{X,r) < 4p(l) r"^ 
A— j-oo 

Proof. Any crossing point is in £{X,r) and so by identity ([2]) 

EN{X,r) <2ttx 27r"^p(A,r) < oo 

and the result follows from (j37|) . ■ 

The following general version can be proved similarly. 

Proposition 19 Lete > and let Ki, K2 be compact sets whose e -neighborhoods 
Ki,K2 are disjoint. For zi £ Ki,Z2 G K2 let TZs{zi,Z2) be the subroute 
of TZ{zi, Z2) crossing from the boundary of Kf to the boundary of Let 
N(X,e,Ki,K2) be the number of distinct paths amongst the sei {7^£(.^i, .^2) ■ 
G E{X)n Ki}. Then 

E lim N{X,e,Ki,K2) < 00. 

A— >oo 

6.2 A bound on normalized length 
Lemma 20 i < 2p(l). 

Proof. Define 

ns{^, e) = m, e) n (disc(e, 6) u disc(e', 6)) 
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in words, the part of the route that is within distance S from one or both 
endpoints. Then define 

'5^(A) = U5,5,g=(;,)7^^(e,^')• 

Note that clearly 

S{X)\Si{X)<Z£{X,l). (38) 

By considering A = 1, 

£<p{l,l) + L{Si{l)) 
where t{-) denotes edge-intensity. Now write 

.(5i(l)) = ^.(521-^(1) \52-.(l)). 
fc>i 

For fixed k > 1, scaling by 2^ gives 

i{S,.^.{l)\S,-.{l)) = 2-'=.(52(2-2^)\5i(2-2'^)) 

< 2-^p{2'^'',l) by §S\). 

So 

^< J;2-v(2-2^l)< j;2-v(i). 

fc>0 fc>0 



6.3 The network £{oo,r) of major roads 

Intuitively, the point of assumption ([20|) and the scaling relation (j37|) is 
that we can define a proces £{oo,r) := UA<oo^(A,r) which must have edge- 
intensity p{r) = p(l)/r, and that in results like Proposition[T8]we can replace 
lim>,__^oo -^(A, r) by N{oo,r). We don't want to give details of a completely 
rigorous treatment, but let us just suppose we can set up £{oo,r) as a 
random element of some suitable measurable space, as we did for 5(A) in 
section 12. 3i 

The conceptual point is that 5(A) and £{X,r) depend on the external 
randomization, that is on the fact that we were studying a SIRSN via the 
random points H(A), but as outlined below £{oo,r) doesn't depend on such 
external randomization. Intuitively this is simply because U;vH(A) is dense 
in M?; we outline a measure-theoretic argument below. 
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Proposition 21 The FDDs (span(2;i, . . . , Zk)) of a SIRSN can be extended 
to a joint distribution, of these FDDs jointly with a random process £*{oo, r), 
such that, for any space-time PPP (H(A),0 < A < oo) independent of the 
FDDs, we have £*{co,r) := [Jx<oo£{^-,''') 

Outline proof. For a suitable formalization of "random subset of M^" we 
have the implication 

if ^1 and A2 are i.i.d. random subsets, and if A\ U A2 —a.s. 
A' = ^1, then Ai = A a.s. for some non-random subset A 

and then the corresponding "conditional" implication 

if Z is a random element of some space, if ^1 and A2 are random 
subsets conditionally i.i.d. given Z, and if ^iU^2 ^a.s. A' where 
{Z,A') = {Z,Ai), then ^1 = ^ a.s. for some Z-measurable 
random subset A. 

So take two independent space- time PPPs H1(A),H2(A) and use a measure- 
preserving bijection [0, 00) U [0, 00) — )• [0, 00) to define another space-time 
PPP H'(A) in terms of and H^. The associated networks satisfy 

£\oo,r)u£^{oo,r) =£'{oo,r) = £\oo,r) 

and this holds jointly with the FDDs of the SIRSN. Since £^{cc,r) and 
iS^(oo,r) are conditionally i.i.d. given the SIRSN. Proposition 1211 follows 
from the general "conditional implication" above. 

6.4 Transit nodes and shortest path algorithms 

Here we make a connection with the "shortest path algorithms" literature 
mentioned in section II. 4[ 

Fix h and take the square grid of lines with inter- line spacing equal to h. 
Define Th to be the set of points of intersection of £{00, h) with that grid. 

Lemma 22 (i) Th has point-intensity 47r~^/i~^p(l). 

(a) For each z G E? there is a subset of Th, of mean size 247r^^p(l), 
and with \z' — z\ < for each z' £ Tz, such that for each pair zi, Z2 with 

\z2 — zi\ > 3h the route TZ{zi,Z2) passes through some point ofT^^ and some 
point ofTz2- 
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Proof. The grid has edge-intensity 2h^^, so from ([2]) the point-intensity of 
Th is 2tt^^ X p[h) X 2h~^ , and (i) foUows from scahng (j37p . 

For any starting point z consider the closest grid intersection {ih,jh). 
Then z is in some square with corner {ih,jh), say the square [(i — l)h, ih] x 
[jh, {j + Let Tz be the set of points of intersection of £'(00, h) with the 
concentric square Sz = [(i — 2)/i, {i+l)h] x [{j — l)h, (j+2)/i]. This square has 
boundary length 12h and so the mean size of equals 27r~^ x p{h) x 12h = 
2A-K~^p{l). By construction 

|/i < \z' — z\ < 2^^'^h for each z' on the boundary ofS'z 

and in particular for each z' £ T^. \z2 — zi\ > 3h then the squares Sz-^ 
and Sz2 do not overlap, and the points z[ and Z2 at which the route crosses 
their boundaries are in Th- ■ 

Informal algorithmic implications One cannot rigorously relate our 
"continuum" setup to discrete algorithms, but in talks we present the fol- 
lowing informal calculation. For the real-world road network in a country 
we have empirical statistics 

• A: area of country 

• ij: average number of road segments per unit area 

• M = TjA: total number of road segments in country 

• p{r): "length per unit area" of the subnetwork consisting of segments 
on routes with start/destination each at distance > r from the seg- 
ment. 

For a real-world network there is an inconsistency between scale-invariance 
and having a finite number r] of road segments per unit area, but let us imag- 
ine approximate scale-invariance over scales of say 2 - 100 miles, and modify 
a scale-invariant model by deleting road segments of very short length. In 
what follows it is helpful to imagine the unit of length to be (say) 20 miles. 

Fix r. Lemma [22] (with h = r) suggests that in the real- world network 
we can find transit nodes such that there are 0{p{l)) transit nodes within 
distance 0(r) of a typical point. If so then we can analyze the algorithmic 
procedure outlined in section [L4l The local search involves a region of radius 
r and hence with 0{rjr'^) edges. Regarding the time-cost of a single Dijkstra 
search as ci x (number of edges), the time-cost of finding the route to each 
local transit node is O (ci (r/r^)p(l)). Transit nodes have point-intensity 
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0(p(l)/r^), so the total number is 0{Ap{l)/r'^). Regard the space-cost of 
storing a k x k matrix of inter-transit-node routes as C2k'^] so this space- 
cost is O (c2 (p(l)j4/r^)^) . Summing the two costs and optimizing over 
r, the optimal cost is 0(^/^4^^7/2/3^2/3^4/3(1)) ^ 0{cf^cl^^p'^/^{l)M^/^) 
and this ©(M^/^) scaling represents the improvement over the 0{M) scal- 
ing for Dijkstra. The corresponding optimal number of transit nodes is 
0((ci/c2)^/^p2/3(x)M"'^/^). The latter has a more interpretable formula- 
tion. If the only alternative algorithms were a Dijkstra search of cost 
ci X (number of edges) or table look-up of cost C2 x (number of edges) ^ , 
then there would be some critical number of edges at which one should 
switch between them, and this is just the solution m^^i^ of cim„it = C2m'^^if 
So the optimal number of transit nodes is 0{mll-^p'^/^{l)M^^^). 

6.5 Number of singly infinite geodesies 

Write 5* (A) for the spanning subnetwork obtained from 5(A) by adding a 
city at the origin 0. This process inherits the scaling-invariance property 
()12p of 5(A). We know from (|29p that at least one singly infinite geodesic 
from exists. The set of all singly infinite geodesies in 5* (A) from forms 
a priori a tree, because two geodesies that branch cannot re-join, by route 
compatability property (iv) from section [2^21 So consider 

^(A, r) := E(number of distinct points at which some singly infinite geodesic 
in S*{X) from first crosses the circle of radius r). 

What we know in general is 

1 < 9(A,r) < oo; r — )• q{X,r) is increasing; A — )• q{X,r) is increasing 

and the scaling property gives 

<7(A,r) = g(l,rAi/2). (39) 

So the A — )■ 00 limit (7(00, r) := lim;j_j.oo 9(A, r) exists (maybe infinite), and 
the scaling property implies 

(7(00, r) = 0^(00, 1) G [1, 00], < r < 00. 

So consider the property 

0^(00, 1) < 00. (40) 
By applying Proposition [18] with r ~ 1 we see 

<?(oo,l) <4p(l). (41) 
So (I36p implies (I40p . So we have shown the following. 
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Corollary 23 A — t- oo the number of singly infinite geodesies in S*{X) 
from increases to a finite limit number (perhaps a random number with 
finite mean) G. Moreover, if G > 1 then these geodesies branch at 0. 

7 Unique singly-infinite geodesies and continuity 

For a SIRSN, let us call the property G = 1 a.s. (in the notation of Corollary 
[23] above) the unique singly-infinite geodesies property. It is conceivable 
that this property always holds - we record this later in Open Problem [32j 
Uniqueness of geodesies is closely related to continuity of routes TZ{zi,Z2) 
as (21,22) vary, as will be seen in section [721 

7.1 Equivalent properties 

Here we show that several properties, the simplest being (|42p . are equivalent 
to the unique singly-infinite geodesies property. We will give definitions and 
proofs as we proceed, and then summarize as Proposition 1241 

Consider two independent uniform random points Ui,U2 in disc(0, 1). 
By the route-compatability property, the intersection of 7^(0, Ui) and 7^(0, U2) 
is a sub-route from to some branchpoint -61,2 j where either Bi^2 7^ or the 
intersection consists of the single point (in which case, set Bi^2 = 0)- So 
we can define a property 

P(5i 2 = 0) = 0. (42) 

Unique singly-infinite geodesies imply (I42p . Suppose (142 p fails. Then 
there exists e > such that, for independent random points Ui,U2 in 
disc(0,l) \ disc(0,e), their branchpoint satisfies P(-B]^2 = 0) > e. 

Scaling by e~"^,m > 1 and using scale-invariance, for independent ran- 
dom points U^,U^ in disc(0,e~™) \ disc(0, e^~™), their branchpoint 
satisfies P(-Bf2 = 0) ^ ^- It follows that, with probability > e — o(l) as 
m — > 00, there exists points CT^^T of ^(1) (disc(0, e"™) \ disc(0, e-*^^™)) 
such that 

routes n{0,CT) and 71(0, C^) branch at 0. 

So on an event of probability > e this property holds for infinitely many m. 
Then on that event we have G > 1, by compactness within the spanning 
subnetwork 5* (1) (Lemma[T]). ■ 

Next consider the spanning subnetwork 5* (A) on points H(A) U {0}. The 
intersection of all routes 7^(0,^), ^ G ^(A) is a sub-route from to some 
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branchpoint -B(A). So we can define a property 

P(-B(l) = 0) = 0. (43) 
Clearly ([13|) implies (02]); we need to argue the converse. 

(j42p implies (j43p . Suppose (|42p . For each r < c« the intersection of 
routes 7^(0,^), ^ G n disc(0,r) is a subroute vr(l,r) from to some 
branchpoint B{1, r), and by (j42p . scaling and the finiteness of H(l)ndisc(0, r) 
we have 

F{B{l,r) = 0) = 0, each r < oo. (44) 

As r increases the subroute 7r(l, r) can only shrink, and the quantity in ()43p 
is the limit B{1) = lim^^oo -8(1, r). To prove (gS]) it suffices, by (jMD, to 
prove 

B{l,r) is constant for all large r, a.s. (45) 

We may suppose (otherwise the result is obvious) that for some > 4 
the subroute 7r(l,ro) stays within disc(0, 1). As r increases, the only way 
that B{l,r) can change at r is if there is a point ^ G ^(1) H circle(0,r) for 
which the route 7^(0,^) diverges from the existing subroute 7r(l,r— ) before 
the existing branchpoint B{l,r—). If this happens, at ri say, then consider 
the subroute 9{ri) = 7^(0,^) n (disc(0,4) \ disc(0,l)) which has length at 
least 3. Now suppose B{l,r) again changes at some larger value r2. Then 
the corresponding subroute 9{r2) must be disjoint from 0{ri), by the route- 
compatability property. Now the "finite length in bounded regions" property 
([U implies that B{l,r) can change at only finitely many large values of r, 
estabhshing (ji5|) . ■ 

Now make a slight re-definition of B{X), by considering only points outside 
the unit disc. That is, the intersection of all routes 7^(0,^), ^ G \ 
disc(0, 1) is a sub-route 7r(A) from to some branchpoint Bi{X). Using 
scale-invariance it is easy to check that (I43p is equivalent to 

P(5i(A) = 0) = for each A < oo. (46) 

As A increases, the sub-routes 7r(A) can only shrink, and the intersection of 
these subroutes over all A < oo is again a subroute from to some point 
Bi{oo). So we can define a property 

P(-Bi(oo) = 0) = 0. (47) 

Clearly (1471) implies (I46p : we need to argue the converse. 
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(j46p implies (I47p . Suppose (jH]). To prove ([T7|) we essentially repeat the 
argument above, but use assumption (pOj) instead of ([5]). It is enough to show 
that, as A increases, -Bi(A) can change at only finitely many large values of 
A. And we may suppose that for large A the subroute 7f(A) stays within 
disc(0,l/4). If -Bi(A) changes at Ai then there is a point ^ appearing at 
"time" Ai for which 7^(0,^) diverges from the existing subroute 7r(Ai— ) and 
so must cross circle(0, 5/8) at some point z{\i) G £{Xi,3/8) C £^(oo,3/8). 
By route-compatability the points z(Aj) corresponding to different values Aj 
where -Bi(A) changes must be distinct, and then ([^n]) implies i5(oo,3/8) n 
circle(0, 5/8) is an a.s. finite set of points. ■ 

Clearly (|47p implies unique singly-infinite geodesies, by the final assertion of 
Corollary [23l We have now shown a cycle of equivalences. Finally, by scaling 
(j47p is equivalent to the following property, where the notation is chosen to 
be consistent with notation in the next section. Define Q{X,0,B) to be the 
probability that the routes 7^(0,^') to all points ^' E 5(A) n (disc(0, B)y do 
not all first exit disc(0, 1) at the same point. Then (I47p is equivalent to 

lim lim Q{X,0,B) = 0. (48) 

Bt-oo A->-oo 

To summarize: 

Proposition 24 Properties (4^ 1, [43\ ), (44\ >! 23' '^"'^ 17^ are each equiv- 
alent to the unique singly-infinite geodesies property. 

7.2 Continuity properties 

In the previous section we studied properties of long routes from a single 
point. We now consider long routes from nearby points, and in this context 
it seems harder to understand whether different properties are equivalent. 
Suppose, for this discussion, the unique singly-infinite geodesies property 
holds. Then the geodesies from and from 1 = (1, 0) G are either disjoint 
or coalesce; we do not know (Open Problem [32]) whether the property 

the geodesies from and from 1 coalesce a.s (49) 

always holds or is stronger. There are several equivalent ways of saying 
()49p - see the end of this section - but what's relevant now is that it is 
equivalent to the property that, for each A, the geodesies from each point of 
5* (A) n disc(0, 1) coincide outside a disc of random radius i?(A) < oo a.s.. 
So we can then ask whether the property 

R{oo) := lim R{X) < oo a.s. 

A— >oo 
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is implied by property (j49p or is stronger. We restate this latter property 
as (ISnD below. 

For < e < 1 < define Q{\,e,B) to be the probability that the 
routes TZiCC) between points ^ S 5(A) H disc(0, e) and points ^' E 5(A) R 
(disc(0, B)Y do not ah first exit disc(0, 1) at the same point. Note Q{X, e, B) 
is monotone increasing at A increases, and decreasing as B increases or e 
decreases. So we can define 

Q{oo,e,B) := lim Q{X,e,B) 

A— >oo 

and then define a property of a SIRSN 

lim Q(oo,e,B) = (50) 

£4,0,Stoo 

where the limit value is unaffected by the order of the double limit. In words, 
(|5Up says that (with high probability) every route from a small neighborhood 
of the origin to any distant point will first cross the unit circle at the same 
place. Property (j50p implies (149 P and implies form (I48p of the unique singly- 
infinite geodesies property, which is the same assertion for routes from the 
origin only. 

The kinds of properties described above relate to questions about conti- 
nuity of the routes TZ{zi, Z2) as zi, 22 vary, and we will give one such relation 
as Lemma [25] below. 

Consider < r/ < 5 < 1/2 and for points ^ G ^(A) n disc(0,7/) and G 
5(A)ndisc(l, ry) with route 7^(^, C') let TZs{^, ?') be the sub-route between the 
first exit from disc(0,5) and the last entrance into disc(l,5). Let ^{\,ri,5) 
be the probability that the sub-routes TZsiCC) for all ^ G ^(A) n disc(0,r/) 
and all ^' G 5(A) n disc(l,r/) are not all an identical sub-route. As above, 
by monotonicity we can define 

^{00, 7], 6) := lim ^{\,7j,5) 

A—)- 00 

and then define a property of a SIRSN 

lim^(oo,ri,6) =0 V(5. (51) 

In words, (I5ip says that (with high probability) all routes from a very small 
neighborhood of the origin to a very small neighborhood of 1 coincide outside 
of larger small neighborhoods. 

Lemma 25 Property / f5Q|) implies property \51\) . 
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Proof. Choose a such that ar] < 1 < a5. Take the definition of scale by 
a, and use scale-invariance to obtain the following. 

The probability that the sub-routes TZasH, CO for all ^ G 5(a^^A)n 
disc(0,aT/) and all C' G 5(a~^A) n disc((a, 0), a?]) are not all an 
identical sub-route equals 'I'(A,7?, 5). 

When this occurs there are two non-identical sub-routes between circle(0, a5) 
and circle((a, 0), a5), which imply two non-identical sub-routes between circle(0, 1) 
and circle((a, 0), 1). For this to happen, either the defining event for Q(a~^A, ar/, a/2), 
or the analogous event with reference to (a, 0) instead of 0, must occur; oth- 
erwise all routes in question pass through the same points on circle(0, 1) and 
circle((a, 0), 1), contradicting the route-compatability properties of section 
[221 So 

^'(A,?7,(5) < 2Q{a-'^\ar],a/2). 

Letting A — >• oo 

^'(oo,r/, (5) < 2Q{oo,ar],a/2). 

Choosing a = ri~^^'^ establishes the lemma. ■ 

Remark. Lemma [25] is almost enough to prove that, under condition 
(|50p . we have the continuity property 

if (z^ z^) ^ {zi,Z2) then 7^(z]', z^) 7^(zl, za) a.s. (52) 

where convergence of paths is in the sense of section 12.31 To deduce (I52p 
one would need also to show that the lengths of Tl{z^,Z2) H (disc(2;i, e„) U 
disc{z2,£n)) tend to a.s. for all e„ — t- 0. This is loosely related to Open 
Problem [M 

Another property equivalent to (I49|) . Because geodesies either co- 
lalesce or are disjoint, for any countable set of initial points there is some 
set of "geodesic ends" , where each such "end" corresponds to a tree of co- 
alescing geodesies from originating "leaves". By a small modification of 
the proof of Corollary [23l the mean number of such ends from the points 
H(A) n disc(0, 1) is at most 4p(l), so we can let A — )• oo and deduce that 
the number G* > 1 of ends from initial points H(oo) n disc(0, 1) satisfies 
EG* < 4p(l). Then by scale-invariance, for each < r < oo the number of 
ends from initial points H(oo) n disc(0, r) equals G* . So the property 

G* = 1 a.s. 
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is clearly equivalent to property ()49p (plus the unique singly-infinite geodesies 
property). Note that if G* > 1 then there are a finite number of different 
"geodesic trees" each of whose leaf-sets is dense in R'^ - behavior hard to 
visualize. 

7.3 The binary hierarchy model 

Proposition 26 The binary hierarchy model has property i50\) . 

Proof. Consider the last stages of construction of the model in section 
13.71 Rotation and scaling do not affect the property of interest, so it will 
suffice to prove the property in the model 7^t-i ■ Consider the argument 
from "proof of Proposition (9]' in section 13.61 but with large rescalings of B 
instead of small rescalings. Combining this argument with the construction 
of 7^t-i at the start of section 13.71 one can show (details omitted) that the 
set 

A' := {z €z M."^ : z in only finitely many B^} 

has area zero; here B^ := a2iB is the "large" rescaling of the union B := 

Ug^g of the translates of the small subsquare S of the basic 2'^+^ x 
2/1+1 

square G in the Figure 6 configuration. By translation-invariance, this 
implies that a.s. A' . For such a realization there is a random infinite 
sequence i(j) with G a2tU)B, and any singly-infinite geodesic from must 
pass through the corresponding infinite sequence 6j(j) of points determined 
by Figure 6. This establishes the unique singly-infinite geodesic property. 
Moreover lies in some translated square of side e2*'^-'^ and for any 
other point in that square its geodesic must coalesce with the geodesic from 
at or before h^y It is easy to check that the squares ^i(j) eventually cover 
any fixed disc, and this establishes property (|50l) . ■ 

8 Open problems and final discussion 
8.1 Other specific models? 

A major challenge is finding other explicit examples of SIRSN models. Let 
us pose the vague problems 

Open Problem 27 Give a construction of a SIRSN which is "mathemat- 
ically natural" in some sense, e.g. in the sense that there is an explicit 
formula for the distribution of subnetworks span(2;i, . . . , z^). 
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Open Problem 28 Give a construction of a SIRSN which is "visually re- 
alistic" in the sense of not looking very different from a real-world road 
network. 

8.2 Quantitative bounds on statistics 

In designing a finite road network there is an obvious tradeoff between total 
length and the network's effectiveness in providing short routes, so in our 
context there is a tradeoff between i and A := EDi. More generally 

Open Problem 29 What can we say about the set of possible values, over 
all SIRSNs, of the triple (A = EDi,l,p{l)) of statistics of a SIRSN? 

This is a sensible question because each statistic is dimensionless, that is not 
dependent on choice of unit of length - a non-dimensionless statistic would 
take all values in (0, oo) by scaling. 

We have given three results relating to this problem. Proposition 1161 gave 
a crude lower bound on the function ^*(A) defined as the infimum value of 
i over all SIRSNs with the given value of A. 

Open Problem 30 (i) Give quantitative estimates of the function £^{A), 
improving Proposition \16l 

(a) Do "optimal" networks attaining the infimum exist, and (if so) can we 
say something about the structure of the associated optimal networks? 

One might make the (vague) conjecture that for some value of A the optimal 
network exploits 4-fold symmetry in some way analogous to our section [3] 
model, and that for some other value it exploits 6-fold symmetry. 

In section 15.51 we showed ()35p that the overall minimum normalized 
length l^, := infA^*(A) satisfies £^, > \/l/8. The third result was Lemma 
[201 showng e < 2p{l). 

8.3 Traffic intensity 

As mentioned in section [2^31 the conceptual point of i5(oo, r) is to capture the 
idea of the major road - minor road spectrum, and the particular definition 
of £"(00, r) is mathematically convenient because of the scaling property ([37|) 
of the edge-intensity p{r). But from a real- world perspective it seems more 
natural to use some notion of traffic intensity. Given any measure ■0 on 
source-destination pairs {zi,Z2), then length measure Lebi along the routes 
Tl{zi,Z2) in a SIRSN induces a "traffic intensity" measure ip on Ur£{(x>,r). 
The natural measures ip to consider are specified by 
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(i) zi has Lebesgue measure Leb2 on 

(ii) given zi, the measure on z := Z2 — zi has density \z\~^ . 
The action of dc on gives the measure specified by 

(i) zi has measure c~^Leb2 on 

(ii) z := Z2 — z\ has density c^~'^\z\~^ 

(iii) the measure along paths is c~^Lebi 
and this is the measure i^) scaled by c^~^ . 

To make a rigorous treatment, the issue is to show that ■i/' is a locally finite 
measure on £"(00, 1). Heuristically one needs /3 > 2 so that the contribution 
from large \z2 — zi\ is finite, and /3 < 4 so that the contribution from small 
\z2 — zi\ is finite. 

Open Problem 31 Show that, perhaps under regularity assumptions on 
the SIRSN, for 2 < /3 < 4 the construction above gives a locally finite mea- 
sure ip on £{00, 1) and hence on Ur£{oo,r). 

8.4 Technical questions raised by results 

8.4.1 Implications between different properties of a SIRSN 

We have given various results of the form "one property of a SIRSN implies 
another" for which we conjecture the converse is false. In particular, we 
expect there are counter-examples to most of the following, though of course 
this requires constructing other examples of SIRSNs. 

Open Problem 32 Prove, or give a counter-example to: 

(i) l[TS\) implies (2^1 

(ii) the unique singly -infinite geodesies property implies 

(iii) (^Pp implies |5Q|) 

(iv) |22P implies [5D\) . 

8.4.2 Understanding the structure of £{00, 1). 

We envisage £^(00, 1) as looking somewhat like a real- world network of major 
roads, but it is not clear what aspects of real networks appear automatically 
in our SIRSN model. For instance, a priori £{00, 1) need not be connected 
(it might contain a short segment in the middle of a route between two 
points at distance 2 + e apart) but it must contain an unbounded connected 
component (most of a singly-infinite geodesic). 

Open Problem 33 Does £{00, 1) have a.s. only a single unbounded con- 
nected component? 
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8.4.3 Questions about lengths 

Even though we started the whole topic of SIRSNs by considering route- 
lengths, they have played a rather small role in our results, and many ques- 
tions about route-lengths could be asked. 

Open Problem 34 Under what extra assumptions (if any) is it true that, 
for Ui,U2, ■ ■ ■ independent uniform on disc(0, 1), 

Esuplen[7^(0,C/^)] < oo? 

i>l 

The following (intuitively obvious) claim seems curiously hard to prove; 
the difficulty lies in showing that the spanning subnetwork does not have 
(necessarily with low probability) huge length a long way away from the 
square. 

Open Problem 35 Take k uniform random points Zi, . . . , in a square 
of area k and consider the length len[span(Zi, . . . , Z^)] of the spanning sub- 
network random network span(Zi, . . . , Z^). Prove 

E len[span(Zi, . . . , Zk)] ^ Ik as k ^ oo. 
8.5 Alternative starting points for a setup 

We started the whole modeling process by assuming we are given routes 
between points, but one can imagine two different starting points. The first 
involves starting with a network of major roads and then adding successively 
more minor roads, so eventually the road network is dense in the plane. In 
other words, base a model on some explicit construction as r decreases of 
some process (<?(r), oo > r > 0) of "roads of size > r" (in our setup this 
is achieved implicitly by the networks £{oo,r)). Of course this corresponds 
to what we see when zooming in on an online map of the real-world road 
network; the maps are designed to show only the relatively major roads 
within the window, and hence to show progressively more minor roads as 
one zooms in. In talks we show such zooms along with the online "zooming 
in" demonstration j24j of Brownian scaling to illustrate the concept of scale- 
invar iance. 

The second, mathematically abstract, approach is to start with a random 
metric d{z,z') on the plane, and define routes as geodesies. 

But a technical difficulty with both of these approaches is that there 
seems no simple way to guarantee unique routes between a. a. pairs of points 
in the plane - in general one needs to add an assumption of uniqueness. The 
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explicit models constructed in section [3] and outlined in section [J] do use the 
"random metric" idea, but the hard part of the construction is proving the 
uniqueness of routes, even in these simplest models we can imagine. It is 
perhaps remarkable that our approach, taking routes as given with only the 
route-compatability property but with no explicit requirement that routes 
be minimum-cost in some sense, does lead to some non-obvious results. 

8.6 Empirical evidence of scale-invariance? 

For real-world road networks, can scale-invariance be even roughly true over 
some range of distance? We mentioned one explicit piece of evidence (or- 
dered segment lengths) in section [T3| there is also evidence that mean route 
length is indeed roughly proportional to distance, though this is also consis- 
tent with other (non scale-invariant) models [6]. 

An interesting project would be to study the spanning subnetworks on 
(say) 4 real-world addresses, whose positions form roughly a square, ran- 
domly positioned, and find the empirical frequencies with which the various 
topologically different networks appear. Scale-invariance predicts these fre- 
quencies should not vary with the side-length of square; is this true? 

8.7 Other related literature 

8.7.1 Hop count in spatial networks 

There has been study of spatial networks with respect to the trade-off be- 
tween total network length and average graph distance (hop count), instead 
of route-length. See |23j for a recent literature survey and empirical analysis. 

8.7.2 Continuum random trees in the plane 

Existence of continuum limits of discrete models of random trees has been 
conjectured, and studied non-rigorously in statistical physics, for a long time, 
and since 2000 spectacular progress has been made on rigorous proofs. For 
three models of random trees (uniform random spanning tree on Z^, minimal 
spanning tree on (with random edge lengths) , and the Euclidean minimal 
spanning tree on Poisson points) , [2J established a rigorous "tightness" result 
and gave sample properties of subsequential limits. A subsequent deep result 
|18j established the existence of a continuum limit in the first model. In these 
limits the paths have Hausdorff dimension greater than 1 so Di = oo a.s.. 
There should be a simple proof of the following, because our definition of 
SIRSN requires KDi = oo. 
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Open Problem 36 In a SIRSN, the subnetwork S {1) cannot be a tree (with 
Steiner points). 

8.7.3 Geodesies in first-passage percolation 

Geodesies in particular models of first-passage percolation have been studied 
in ^19j. It is unclear whether there is any substantial connection between 
the behavior of geodesis in that setting and in our setting. 

8.7.4 A Monge-Kantorovitch approach 

A completely different approach to continuum networks, starting from Monge- 
Kantorovitch optimal transport theory, is developed in the monograph by 
Buttazzo et al. Their model assumes 

(i) some continuous distribution of sources and sinks 

(ii) an a priori arbitrary set S representing location of roads 

(iii) two different costs-per-unit-length for travel inside [resp. outside] S. 
An optimal network in one that minimizes total transportation cost for a 
given cost functional on S. It is shown that, under regularity conditions, 
the optimal network is covered by a finite number of Lipschitz curves of 
uniformly bounded length, although it may have even uncountably many 
connected components. But this theory does not seem to address statistics 
analogous to our A and £ in any quantitative way. 

8.7.5 The method of exchangeable substructures 

The general methodology of studying complicated random structures by 
studying induced substructures on random points has many applications [7j. 
In particular, the Brownian continuum random tree \3} provides an analogy 
for what we would like to see (Open Problem \27\i in some "mathematically 
natural" SIRSN - see e.g. the formula (13) therein for the distribution of 
the induced subtree on random points - though that is in the "mean-field" 
setting without any d-dimensional geometry. 

8.7.6 Urban road networks. 

There is scattered literature on models for urban road networks, mostly with 
a rather different focus, though [IT] has some conceptual similarities with 
our work. 



52 



8.7.7 Dynamic random graphs. 

Conceptually, what we are doing with routes TZ{zi,Z2) and subnetworks 
5(A) is exploring a given network. This is conceptually distinct from using 
sequential constructions of a network, a topic often called dynamic random 
graphs [13j, even though the particular "dynamic Gabriel" model outlined 
in section 18.11 does fit the "dynamic" category. 

Acknowldgements My thanks to Justin Salez for details of the proof of 
Lemma [TT| to Wilfrid Kendall for ongoing collaboration, and to Cliff Stein 
for references to the algorithmic literature. 
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A Appendix: A topology on the space of feasible 
subnetworks. 

We first define convergence of routes. Recall a feasible route r(z,z') involves 
line segments between points (zj) which we will call turn points of the route. 
Given e < \z' — z\/2 there is (starting from z) a last turn point z^^-j before 
the route r(z, z') first exits disc(z, e) and there is (starting from z') a last 
turn point before the reverse route r{z',z) first exits disc(z',e). Define 

r{z{n),z'{n)) — ^ r{z,z') 

to mean 

(i) z{n) — )• z, z'{n) z' ^ z. 

(ii) For each e < \z' — z\/2 such that circle(z,e) and circle(z',e) do not 
contain any turn point of r(z,z'), writing the turn points of the subroutes 
''(%)W>2(£)("')) as (7/o(n),yi(n), . . . ,yfc(n)), we have 

(yo(n),yi(n), . . . ,yfc(n)) {yo,yi, . . . ,yk) 

the limit being the turn points of the subroute r(2(£-), z'^i,^), where k is finite 
from the definition of feasible route. 

(iii) The total lengths L(^^^{n) of r{z{n), z' {n))n{disc{z, e)Lldisc{z' , e)) satisfy 

lim limsupL(£)(n) = 0. 
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Despite its inelegant formulation, this seems the "natural" notion of conver- 
gence. 

Now we specify, in a way analogous to (i-iii) above, what it means for 
a sequence s(n) of feasible subnetworks on locally finite sets z(n) = {z^{n)} 
to converge to a limit subnetwork s on z. 

(i) We need z(n) to converge to z in the usual sense of convergence of simple 
point processes [12j. This is equivalent to saying that if we take any R such 
that circle(0, R) contains no point of z, then we can label the points of z(n)n 
disc(0,i?) as {z^ (n) , . . . , (n)) in such a way that {z^ (n) , . . . , z^ (n)) — )■ 
(z^, . . . , z^) , the limit (here and in analogous assertions below) being the 
points of z n disc(0, R). 

(ii) Take R and (z^, . . . , z^) as above and take e < ^ mini<j<j<x \ — 
such that Ui<j<ii-circle(2;*, e) does not contain any turn point within s. Then 
we can label the turn points of Ui<j<i^r(2^^-j(n), z;^^^(n)) as (y"(n), I < u < 
L) in such a way that 

(?/"(n), l<u<L)^ (y", 1 < ti < L) 

{y^ (n) , y'" (n)) is an edge-segment of route r(z*(n), z-'(n)) iff {y'^,y'") is an 
edge-segment of route r{z'^,z^). 

(iii) For each 1 < i < j < K the routes r{z^[ri), z^ [n)) satisfy (iii) above. 

(iv) Lemma [U (i) implies that given R, the following quantity (referring to 
the subnetwork s) is finite: 

R* := min{r : |J r{zi, Zj) n disc(0, R) =3(1 disc(0, R) } 

Zi,Zj ^disc{0,r) 

(that is, each edge of s within disc(0, R) is part of some route between 
endpoints in disc(0, i?*)). We require 

lim sup i?*(n) < oo for each R < oo. 

n 

We have described sequential convergence within the space of feasible 
subnetworks. It is routine to show this is convergence in some complete 
separable metric space, but we won't pursue such theory here. 
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